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Variations of the Fine Structure Constant in 

Space and Time 



David Fonseca Mota 
Summctry 

This thesis describes a detailed investigation of the effects of matter inhomo- 
geneities on the cosmological evolution of the fine structure constant. 

In chapter 1, we briefly describe the observational and theoretical motivations 
to this work, and we review the standard cosmological model. We also review the 
Bekenstein-Sandvik-Barrow-Magueijo (BSBM) theory for a varying fine structure 
constant, a. 

Assuming a Priedmann universe which evolves with a power-law scale factor, 
a = t", in chapter 2, we analyse the phase space of the system of equations that 
describes a time- varying a, in a homogeneous and isotropic background universe. We 
classify all the possible behaviours of a{t) in ever-expanding universes with different 
n and find exact solutions for a{t). In general, a will be a non-decreasing function 
of time that increases logarithmically in time during a period when the expansion 
is dust-dominated (n = 2/3), but becomes constant when n > 2/3. a tends rapidly 
to a constant when the expansion scale factor increases exponentially A general set 
of conditions is established for a to become asymptotically constant at late times in 
an expanding universe. 

In chapter 3, using a gauge-invariant formalism, we derive and solve the linearly 
perturbed Einstein cosmological equations for the BSBM theory. We calculate the 
time evolution of inhomogeneous perturbations of the fine structure constant, ^ on 
small and large scales with respect to the Hubble radius. In a radiation-dominated 
universe small inhomogeneities in a decay on large scales but on scales smaller than 
the Hubble radius they undergo stable oscillations. In a dust-dominated universe, 
small inhomogeneous perturbations in a approach a constant on large scales and on 
small scales they grow as t^^^, and ^ tracks ^ . If the expansion accelerates, as in 
the case of a A or quintessence-dominated phase, inhomogeneities in a decay on both 
large and small scales. The amplitude of perturbations in a are much smaller than 
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that of matter or radiation perturbations. We also present a numerical study of the 
non-linear evolution of spherical inhomogeneities in radiation and dust universes by 
means of a comparison between the evolution of flat and closed Priedmann models 
with time-varying a.Various limitations of these simple models are also discussed. 

Chapter 4 is dedicated to the effects of non-linear structure formation in the 
evolution of a. We study the space-time evolution of the fine structure constant 
inside evolving spherical overdensities in a A Cold Dark Matter (CDM) Priedmann 
universe, using the spherical infall model. We show that its value inside virialised 
regions will be significantly larger than in the low-density background universe. The 
consideration of the inhomogeneous evolution of the universe is therefore essential 
for a correct comparison of extragalactic and solar system limits on, and observa- 
tions of, possible time variation of a and other constants. Time variation of a in the 
cosmological background can give rise to no locally observable variations inside viri- 
alised overdensities like the one in which we live, explaining the discrepancy between 
astrophysical and gcochcmical observations. 

In chapter 5, using the BSBM varying-alpha theory, and the spherical collapse 
model for cosmological structure formation, we study the effects of the dark-energy 
equation of state and the coupling of a to the matter fields on the space and time 
evolution of a. We compare its evolution inside virialised overdensities with that in 
the cosmological background, using the standard (A = 0) CDM model of structure 
formation and the dark-energy modification, wCDM. We find that, independently 
of the model of structure formation one considers, there is always a difference be- 
tween the value of alpha in an overdensity and in the background. In a SCDM 
model, this difference is the same, independent of the virialisation redshift of the 
overdense region. In the case of a wCDM model, especially at low redshifts, the 
difference depends on the time when virialisation occurs and the equation of state of 
the dark energy. At high redshifts, when the wCDM model becomes asymptotically 
equivalent to the SCDM one, the difference is constant. At low redshifts, when dark 
energy starts to dominate the cosmological expansion, the difference between a in 
a cluster and in the background grows. 

The last chapter contains a summary of the results obtained in the thesis and a 
discussion of some open problems that require further investigation. 
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Introduction 



. . . contudo fazemos da nossa angustia, 
desta quietude irriquieta que nos atormenta a alma, 
as mil e uma perguntas sem resposta. . . 

- Ernesto Mota - 



1.1 Constants of Nature 

Once upon a time, humankind started to describe the universe we hve in using 
mathematics. When a mathematical formahsm is proposed in order to describe the 
different phenomena that occur in our universe, a theory in physics is set. Any 
theory which claims to explain a set of events in nature has to be confronted with 
the empirical observations. Hence, having the formal structure of a given physical 
theory, there is the need to quantify it in order to test it. The quantification of 
a given physical model consist in the introduction of numbers. But this raises a 
question: How shall we choose these numbers and to which quantities shall we 
attribute them? Theories in physics have several free parameters of which value 
cannot be calculated or deduced from its formal structure. These free parameters 
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are the fundamental constants of the theory. We need to attribute a value to these 
quantities, and the value is chosen in order to match the quantitative predictions 
of the theory with the observations and experiments performed. Why are these 
quantities fundamental? Because their values cannot be calculated, only measured. 
Why constant? Because every experiment yields the same value independently of 
time, position, temperature, pressure, etc. 

But not all the constants in physics are fundamental. If one conducts an experi- 
ment with greater precision and/or in an environment more extreme than previously 
considered, one might observe variations of the constants. This may precipitate 
physicists to search for a more basic theory which explains the origin and the value 
of the (then non- fundamental and variable) constant. The status of a constant de- 
pends then on the considered theory and on the observer measuring them, that is, 
on whether this observer belongs to the world of low-energy quasi-particles or to the 
high energies one. For example, the density of a given material was once regarded 
as constant, but was later found by experimenters to vary with temperature and 
pressure. 

1.1.1 Why study variations of the fundamental constants? 

By definition, the constants of nature are quantities which have no dependence on 
other measured parameters, it is not surprising that despite the enormous advances 
in particle physics we still have no idea why the constants of nature take the values 
they do. This might be something we will never be able to answer, but we have 
learnt from the past, that what we call today a constant might not be so in the 
future. 

The hypothesis of the constancy of the fundamental constants plays an important 
role in astronomy and cosmology where the redshift measures the look-back time. 
Ignoring the possibility of varying constants could lead to a distorted view of our 
universe and if such variations are established corrections should be applied. It 
is thus important to investigate that possibility, especially as the measurements 
become more precise. 

The study of the variations of the constants of nature also offers a new link 
between astrophysics, cosmology and high-energy physics complementary to early 
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universe cosmology. In particular, the observation of the variability of the funda- 
mental constants constitutes one of the few ways to test directly the existence of 
extra-dimensions and to test high energy physics models. 

From all the possible fundamental constants of nature we will focus on the so 
called fine structure constant, a. The fine structure constant is also regarded as the 
coupling constant of the electromagnetic interactions. The fine structure constant 
can be derived from other constants as follows [6]: 

a = - (1.1.1) 

where c is the speed of hght in vacuum, h = h/27r is the reduced Planck constant, e 
is the electron charge magnitude and eo is the permitivity of free space. The value 
of a measured today on earth is ao 1/137.035 [6]. 

In this thesis we will study the theoretical possibilities of time and spatial vari- 
ations of the fine structure constant along the history of the universe. 

Time variations of a can be measured using the 'time shift density parameter': 

Aa ^ a{z) - ao ^ 2) 

a ao 

where a{z) is the value of the fine structure constant at a redshift z. Similarly, 
spatial variations of a can be studied using the 'spatial shift density parameter': 

5a ^ a^j^{z) - a^^{z) ^ 
a oij:^(z) 

where a^j^ j is the value of the fine structure constant in a two different regions in 
the universe, for instance xi and X2 can represent two different galaxies. 

But why shall we look at possible variations of the fine structure constant? The 
next two sections describe both observational and theoretical motivations for this 
work. 
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1.2 Observations and Constraints on Variations 
of the Fine Structure Constant 

Many methods were used to observe possible variations of the fine structure constant 
[7]. A general conclusion of those experiments and observations, is that there is no 
variation of a, or at least the results are consistent with the hypothesis of a null 
variation. If there are any time or spatial variations of a, these will be very small. 
Hence, up to nowadays, we basically have upper limits that constraint variations 
of the fine structure constant. Nevertheless, there is at least one observation which 
explicitly gives a non-null result for variations of a. As we will see, this observation 
has a very unique feature: It has observed variations in the fine structure constant 
in a redshift range of ^ = 0.5 — 3.5. 

When taking conclusions from experiments and observations on variations (or 
none) of the constants of nature we should always bear in mind that those mea- 
surements are set on a given time-scale. Hence we need to be very careful when 
extrapolating those results into an absolute conclusion. For instance, in the case of 
the fine structure constant, one expect to be able to constrain a relative variation of 
a in a time scale of 10^ years {z ~ 0.3) using geochcmical constraints [8]; 10^ — lO^'' 
years {z ^ 3) using astrophysical (quasars) methods [5]; 10^° years {z ~ 1000) using 
cosmological methods (Cosmic Microwave Background Radiation) [9]; and 1 — 12 
months using laboratory methods [10]. It is then not a straight forward task to 
interpret and take absolute conclusions from most of the observations performed 
until today. For example, if we observe no time variations of a using geochemical 
methods, we can only conclude that the fine structure constant did not vary from a 
redshift z — 0.3 up to today. 

Another common feature to the observational and experimental results is the 
presence of degeneracies. Usually, there are other parameters in the model, used 
to study the physical phenomenon, that may infiuence the results. Due to that, 
the conclusions taken are then based on some assumptions, which usually lead to 
an extra error bar besides the systematic experimental errors. The problem of 
the degeneracies is not an easy one, since such analyses are also dependent on our 
understanding of the fundamental interactions [11, 12]. For instance, grand unifying 
theories predict that all the constants should vary simultaneously. One needs to 
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further study the systematic errors of the observations and to propose and reahse 
new experiments which are less model dependent [13, 14]. 

In this section, we shall not attempt to give a full review, deferring the reader 
to [7] and the references therein. 

1.2.1 Geochemical constraints 

All the geological studies are on time scales of order of the age of the Earth, meaning 
that constraint variations of the fine structure constant in a range of redshifts of 
z ^ 0.1 — 0.15, depending on the values of the cosmological parameters. 

Atomic clocks are one of the principal methods we have to measure possible 
variations of the fine structure constant on Earth. The latest constraint on possible 
variations of a, using atomic clocks, is [15]: 

< 4.2 X 10~^^ yr-^ . (1.2.1) 

Notice that, although the constraint is very tight, it is only valid for a very short 
period of time. Usually a few months, the time the experience lasted. 

One of the oldest observations, with non-null results, looking for variations of the 
fine structure constant came from the analysis of the Oklo phenomenon - a natural 
nuclear fission reactor which operated in Gabon, Africa, at approximately 1.8 billion 
years ago [16, 17, 18, 19]. The first analysis were done by Shlyakhter [20] and an 
upper limit on time variations of a resulted to be: Aa/a < 10~^ over a period of 
1.8 billion years. However that study provided two possible results: 

- = (-0.2 ± 0.8) X 10~^^ per year (1.2.2) 
a 

- = (4.9 ± 0.4) X 10"^^ per year (1.2.3) 
a 

which come from two possible physical branches [16, 20, 21]. Note that only the first 
solution is consistent with a null variation of the fine structure constant, whereas 
the other implies a higher value of a in the past. 

Recently many others analysis were performed using the Oklo phenomenon [8, 
22, 21]. Up to today, the tighter constraint using new samples from the Oklo reactor 



a 
a 
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is: 

— ^ (-0.04 ± 0.15) X 10-^ (1.2.4) 

a 

These limits should be considered with care, since in general there are consider- 
able uncertainties regarding the complicated Oklo environment and any limits on a 
derived from it can be weakened by allowing other interaction strengths and mass 
ratios to vary in time as well [23, 21]. 



Constraints at slightly higher redshifts 

Another geochemical process, which is sensitive to changes in a, is the /3-decay 
rate of ^^^i?e [24]. Updating the work of Dyson [25, 26], Olive et al. [21] claim 
that, improved measurements of the ^^^Re /3— decay rate (derived from meteoritic 
^^"^ Re/^^'^Os abundances ratios), imply |Aq;/q;| < 3 x 10"^ over the 4.65 Gyr history 
of the solar system [z 0.45). However, Uzan [7] notes that the above result assumes 
that variations of the /3-decay rate depend entirely on variations of a and not on 
possible variations of the weak couphng constant, see also [27]. Very recently, Olive 
et al. [28], have rederived this constraint, in models where all gauge and Yukawa 
couplings vary independently, as would be expected in grand unification theories. 
The new constraint is | Aq;/q;| < (—8 ± 8) x 10~^. In spite of the potential problems 
of these experiments, the study of this phenomenon is quite important. It provides 
constraints to variations of a at slightly higher redshift than the usual geochemical 
constraints, and at the same time, it give us a hint on the possible evolution of the 
fine structure constant inside our local system of gravity. 

1.2.2 Astrophysical observations 

Astrophysical observation are probably one of the most important types of experi- 
ments which can be performed in order to study the possibility of variations of the 
fine structure constant in cosmological scales. 

Observing the quasar absorption spectra give us information about the atomic 
levels at their position and time of emissions. It is usually assumed, in those ob- 
servations, that all transitions have the same a dependence. So any variation of 
the fine structure constant will affect all the wave lengths by the same factor. An- 
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other assumption is the independence of this uniform shift of the spectra from the 
Doppler effect due to the motion of the source or to the gravitational field where 
it sits. The idea is to compare different absorption lines from different species to 
extract information about different combinations of the constants at the time of 
emission and then to compare them with the laboratory values (see [29] for a de- 
tailed description). A quite surprising result was obtained from these observations: 
an explicit time variation of the fine structure constant was found. The value of a, 
in a redshift range oi z — 0.5 — 3.5, was smaller than the standard value on earth cto 
[5, 30, 31, 32, 33] (see however [34] for another interpretation of the results). The 
latest value for the time shift parameter of the fine structure constant is [35] : 

= -0.543 ± 0.116 X 10"^ (1.2.5) 

a 

for a redshift range of z — 0.5 — 3.5. 

Many other astrophysical observations were performed [36, 37, 38, 39, 40, 41], but 
none of them has explicitly indicate the existence of variations of the fine structure 
constant, that is, they are just upper bounds. This show us how relevant and 
important are the series of results obtained from the quasar absorption spectra. 
Such a non-zero detection, if confirmed, would bring tremendous implications on 
our understanding of physics, for instance, it would be interesting to investigate if 
such a variation is compatible with the test of the universality of free fall [42] . 

Note that, a feature of the results from the quasar absorption spectra, is that they 
are 'incompatible' with the geochemical observations if the variation of a is linear 
with time. If both the geochemical and the quasar's absorption spectra observations 
are confirmed, then the theoretical models proposed to explain variations of a will 
need to explain the reason for such a discrepancy between astrophysical observations 
and geochemical ones. In particular, the constraint of Aa/a < 3 x 10^^ at redshift 
z = 0.45, coming from the /9-decay experiments, is very difficult to be compatible 
with an explicit variation of the fine structure constant of Aa/a 5 x 10~^ at 
z — 0.5, coming from the quasar's absorption spectra. 
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1.2.3 Early universe constraints 

Cosmic Microwave Background Radiation 

Time or spatial variations of the fine structure constant affect the power-spectrum of 
the Cosmic Microwave Background Radiation (CMBR) anisotropics [43, 44]. This is 
due to the fact that a different value of a would affect the biding energy of hydrogen, 
the Thompson scattering cross section and the recombination rates. A smaller value 
of a would postpone the recombination of electrons and protons, i.e., the last scatter- 
ing would occur at a lower redshift. It would also alter the baryon-to-photon ratio at 
last-scattering, leading to changes in both the amplitudes and positions of features in 
the power spectrum. The strongest current constraints dX z 1000 from the CMBR 
power spectrum are at the l^aja ~ 10^^ level if one considers the uncertainties in, 
and degeneracies with, the usual cosmological parameters (for instance Q^, Qa, etc) 
[9, 45, 46, 47, 48, 49]. However, note there is a critical degeneracy between variations 
of a and the electron mass mg [50, 51], since the relativistic-quantum corrections 
of the electron mass depend on the strength of the electromagnetic interaction and 
consequently on a. This degeneracy dramatically weakens the current constraints. 

Big Bang Nucleosynthesis 

The theory of Big Bang Nucleosynthesis (BBN) is one of the corner-stones of the hot 
Big Bang cosmology, successfully explaining the abundances of the light elements, D, 
^i^e, '^He and Li. Assuming a simple scaling between the value of a and the proton- 
neutron mass difference a limit can be placed on /\a/a [52]. A similar method can 
be used considering simultaneous variations of the weak, strong and electromagnetic 
couphngs on BBN [53]. However, in general, estimates based on BBN abundances 
of ^ife suffer from the crucial uncertainty of electromagnetic contribution to the 
proton-neutron mass difference. Much weaker limits are possible if attention is 
restricted to the nuclear interaction effects on the nucleosynthesis of D, ^He and 
^Li [54, 55]. Given the present observational uncertainties in the light element 
abundances, the most conservative limits are Aa/a = (3 ± 7) x 10~^. 
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1.3 Theoretical Motivations 

In the previous section we have described the empirical motivations to study varia- 
tions of the fine structure constant. In reality, they were not very convincing, since 
most of the results are basically upper bounds to variations of a and the only explicit 
results, that claimed a lower value of a in the past, was presented by a sole group 
of collaborators in a series of experiments [35, 5, 30, 31, 32, 33]. 

The situation is completely different from the theoretical point of view. We 
will see that there are several important reasons why we should seriously take into 
account the possibility of a varying-a. Most of them come from theoretical models 
which try to give a more complete description of the universe. 

These theories have strong observational predictions, for instance, in the case of 
a time- varying a, they allow us to evaluate the effects of a varying a on free fall 
which leads to potentially observable violations of the weak equivalence principle 
[56, 57, 58, 42]. They also allow us to investigate whether or not other cosmological 
observations like the Cosmic Microwave Background Radiation arc consistent with 
the variations of a to fit the quasar observations [46, 9, 51]. 

Once again we shall not attempt to give a full review here, deferring the reader 
to [7] and [59] for a more thorough discussion. 

1.3.1 Multi-dimensional theoretical models 

There are very strong theoretical motivations to study time and spatial variations 
of the fine structure constant. The reasons come essentially from novel theories 
which claim to be more 'fundamental' and so are candidates to be the 'Theory of 
Everything'. The best candidates for unification of the forces of nature, in a single 
unified theory of quantum gravity, only seem to exist in finite form if there are more 
dimensions of space than the common four that we are familiar with. This means 
that the true constants of nature are defined in a higher dimensional world and the 
constants we observe are merely effective values, that is, they are a four-dimensional 
projection of the real fundamental constants. So the coupling constants we observe 
and measure may not be fundamental and so do not need to be constant. 

Since we do not observe any other dimensions than the common four, there is a 
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strong indication that, if those extra-dimensions exist, they are indeed very small. 
A general feature of the multi-dimensional theories, is that the size of the extra- 
dimensions is associated to a scalar field. This field is a dynamical quantity, which 
has the desired property of making the extra-dimensions small and stable (in the 
sense that their size will not change) . The reason why we desire the stability of the 
extra-dimensions, is related to the fact that in most of the cases that scalar field will 
be coupled to the matter fields. Any time or spacial variations of this field will then 
be 'seen' as a variation of the interactions' coupling constants. Up to nowadays, no 
mechanisms of stabilising the scalar field, was found. Hence, until this mechanism 
is to be found, a common feature of the multi-dimensional theories is the existence 
of slow changes in the scale of the extra-dimensions which would be revealed by 
measurable changes in our four- dimensional "constants" . 

Multi-dimensional theories also predict relations between different constants, due 
to relations between the coupling of dilaton type fields and matter, see for instance 
[60, 61, 62, 63] for more complex effects then just the variation of the coupling 
constants. 

One of the earliest attempts to create a single unified theory was presented by 
Kaluza [64] and Klein [65] . They considered a five-dimensional model with the aim 
to unify electromagnetism and gravity (for a review see [66] ) . The fifth dimension is 
compactified in order to turn it small. The compactification results into a coupling 
between the scalar field, which is related to the radius of the fifth dimension, and 
the matter fields. Variations of the scalar field will be reflected as variations of all 
the coupling constants. Various functional forms for monotonic time variations have 
been derived where usually the four-dimensional gauge couplings vary as the inverse 
square of the mean scale of the extra dimension (for instance see [60, 67, 68]). Non- 
monotonic variations of a with time where analysed in [69] and the requirements 
for a self-consistent relations, if there are simultaneous variations of the different 
coupling constants, were discussed. Other five-dimensional theories, based on the so 
called Brane- World models, were also used to study variations of the fine structure 
constant. See for instance [7, 59, 70, 71] and references therein. 

A much more recent and popular candidate to a theory of unification is string 
theory. Superstrings theories offer a theoretical framework to discuss the value of 
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the fundamental constants since they become expectation values of some fields [7]. 
One of the definitive predictions of superstring theories is the existence of a scalar 
field, the dilaton, that couples directly to matter [72] and whose vacuum expectation 
value determines the string coupling constant [73]. The four-dimensional couplings 
are determined in terms of a string scale and various dynamical fields, for instance, 
the dilaton and the volume of compact space. Due to the coupling between the 
dilaton and all the matter fields, we expect the following effects: a scalar mixture 
of a scalar component inducing deviations from general relativity, a variation of the 
coupling constants, and a violation of the weak equivalence principle [7]. Various 
string theories will exhibit different variations of the effective couplings; in partic- 
ular, the cosmological evolution of the fine structure constant will depend on the 
form of the coupling between the dilaton and the matter fields which interact elec- 
tromagnetically. The undefined structure of the coupling gave origin to many low 
energy models and approaches to study time and spatial variations of the funda- 
mental constants and its relation to the extra-dimension radius (see for instance 
[52, 74, 75, 76, 77]). In [58] it was shown that cosmological variations of a may 
proceed at different rates at different points in space-time (see also [78, 79, 80]). 
This mean that variations of the fine structure constant may not only be in time 
but in space as well. 

1.3.2 The Bekenstein model for varying-a 

Independently of any extra-dimensional model, Bekenstein [81, 82, 83] formulated 
a framework to incorporate a varying fine structure constant. Working in units in 
which h and c are constants, he adopted a classical description of the electromagnetic 
field and made a set of assumptions to obtain a reasonable modification of Maxwell 
equations to take into account the effect of the variation of the elementary charge, 
e. His eight postulates are: (i) For a constant a electromagnetism is described by 
Maxwell theory and the coupling of the potential vector to mater is minimal, (ii) 
The variation of a results from dynamics, (iii) The dynamics of electromagnetism 
and a can be obtained from an invariant action, (iv) The action is locally invariant 
(v) The action is time reversal invariant, (vi) Electromagnetism is causal, (vii) The 
shortest length is the Plank length, (viii) Gravitation is described by a metric which 
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satisfies Einstein equations. 

Assuming that the charges of all particles vary in the same way, one can set 
e = eQe{x^) where e{x^) is a dimensionless universal field and Cq is a constant 
denoting the present value of e. This means that some well established assumptions, 
hke charge conservation, must give away [84]. Still, the principles of local gauge 
invariance and causality are maintained, as is the scale invariance of the field e. 

Since e is the electromagnetic coupling, the e field couples to the gauge field as 
e^ju in the Lagrangian and the gauge transformation which leaves the action invari- 
ant is eAfj^ — > eAu + d^x, rather than the usual An — > A^ + d^x- The electromagnetic 
tensor generalises to 

F^,. = e-' (d^ieA,) - d^ieA^)) (1.3.1) 

which reduces to the usual form when e is a constant. The electromagnetic action 
is given by 

Sem = j dx'V^F^^F'^'' (1.3.2) 

and the dynamics of e are controlled by the kinetic term action 

Ihc f . , — d,,ed>^e 



j dx'V^^^ (1.3.3) 



where I is a length scale of the theory, introduced for dimensional reasons , and 
which needs to be small enough to be compatible with the observed scale invariance 
of electromagnetism. This constant length scale gives the scale down to which the 
electric field around a point charge is accurately Coulombic {Ipi < I < 10"^^ — 
10~^^cm) to avoid confiict with experiments [81, 85]. 

Ohve and Pospelov [57] generahsed the Bekenstein model to allow additional 
coupling of a scalar field = -Bf(0) to non-baryonic dark matter (as first proposed 
by Damour [86]) and cosmological constant, arguing that in certain classes of dark 
matter models, and particularly in supersymmetric ones, it is natural to expect that 
(p would couple more strongly to dark matter than to baryons. 

The formalism developed by Bekenstein was also in the braneworld context 
[87, 70] and Magueijo [88] studied the effect of a varying fine structure constant 
on a complex scalar field undergoing an electromagnetic U{1) symmetry breaking 
in this framework (see also [89, 90, 91]). Also, inspired by the Bekenstein model. 
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Armendariz-Picon [92] derived the most general low energy action including a real 
scalar field that is local, invariant under space inversion and time reversal, diffeo- 
morphism invariant and with a C/(l) gauge invariance. 

1.3.3 Other investigations 

When discussing variations of the fundamental constants, careful distinction should 
be made between dimensional and dimensionless constants. The reason is the fact 
that every experimental measurement consist in comparing the value of two quan- 
tities. Hence, any measurement of a dimensional constant needs to be followed by 
the units chosen [24]. This fact has important consequences in theoretical models 
which describe variations of a dimensionless quantity, if the later is obtained via a 
combination of other (dimensional) quantities. For instance, this degeneracy, among 
combination of dimensional constants, lead Dirac [93, 94] to propose the 'Large Num- 
bers Hypothesis'. It sates that the existence of certain large dimensionless numbers 
which arise in combinations of some cosmological numbers and physical constants 
was not a coincidence but a consequence of an underlying relationship between them. 

Speaking of variations of dimensional constants is then somewhat ambiguous, 
since observations and experiments can only measure dimensionless quantities. This 
means that we would not be able to experimentally distinguish variations of the fine 
structure constant resulting form a variation of the speed of light or in the electron 
charge (see however [95]). Bekenstein's model assumed that variations of the fine 
structure constant were due to variations of the electron charge. Another approach is 
to attribute the change in the fine structure constant to a varying light propagation 
speed [96, 97, 98, 99, 100, 101]. The motivation for such theories are their ability to 
solve the standard cosmological problems, usually solve by infiation[102, 103, 104, 
105]. For instance, the horizon problem is trivially solved by claiming a much higher 
light speed in the early universe. 

1.4 Summary 

Observations of quasar absorption-line spectra has been found to be consistent 
with a time variation of the value of the fine structure constant between redshifts 



1.5 Thesis Aim 



z — 0.5 — 3.5 and the present. The entire data set of 128 objects gives spectra consis- 
tent with a shift of a with respect to the present value [35, 5, 31, 32]. Nevertheless, 
extensive analysis has yet to find a selection effect that can explain the sense and 
magnitude of the relativistic line-shifts underpinning these deductions. Adding to 
that, astronomical probes of the constancy of the electron-proton mass ratio have 
reported possible evidence for time variation, [106] , but as yet the statistical signif- 
icance is low. 

Motivated by these observations, there has been considerable theoretical inves- 
tigation of the cosmological consequences of varying a, [16]. In particular, in the 
theoretical predictions of gravity theories which extend general relativity to incor- 
porate space-time variations of the fine structure constant. These have been primary 
formulated as Lagrangian theories with explicit variation of the velocity of light, c, 
[98, 96, 99, 100, 105], or of the charge on the electron, e, [82, 107, 108, 109, 110]. A 
range of variant theories have been investigated with attention to the possible parti- 
cle physics motivations and consequences for systems of grand and partial unification 
[61, 111, 57, 69, 58, 112, 113, 92, 114, 115]. 

Varying-a theories offer the possibility of matching the magnitude and trend 
of the quasar observations and to investigate whether or not other cosmological 
observations are consistent with the small variations of a. They are also of particular 
interest because they predict that violations of the weak equivalence principle [57, 58] 
should be observed at a level that is within about an order of magnitude of existing 
experimental bounds [95, 56, 100]. 

1.5 Thesis Aim 

As we saw in the previous sections, studies on the time and spatial variations of the 
fine structure constant, a, are motivated by two main aspects: 

• The first, are the recent observations of small variations of relativistic atomic 
structure in quasar absorption spectra, which suggest that the fine structure 
constant, was smaller at redshifts z = 0.5—3.5 than the current terrestrial value 
ao = 7.29735308x10-3, with Aa/a = {a(;z)-ao}/ao = -0.543±0.116x 10"^. 

• The second is the fact that the current theoretical models, candidates to 
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a grand unified theory of quantum gravity, require the existence of extra- 
dimensions beyond the common four we are used to. In these models, the 
radius of the extra-dimensions is associated to a scalar field which is coupled 
to the matter fields. Any change in this scalar field refiects variations of the 
coupling constants, in particular the fine structure constant. 

We also saw in the previous sections, that several theories have been proposed to 
investigate the implications of a varying fine structure constant. A common feature 
to all of them is the existence of a field responsible for variations of a. This field is 
coupled to the matter fields, at least, to the ones which interact electromagnetically. 
Due to the coupling between the field and matter, it is then natural to wonder if 
the existence and evolution of inhomogeneities of the later will affect the evolution 
of a. This will be the main objective of this thesis: to investigate how the evolution 
of the inhomogeneities in the matter fields afi^ect the time and spatial evolution of 
the fine structure constant along the history of the universe. 

In order to investigate the effects of the evolution of the matter inhomogeneities, 
on the time and spatial variations of the fine structure constant, we need to first to 
study the evolution of a homogeneous and isotropic universe and of its components. 
In the next section we briefiy review the the Standard Cosmological Model and the 
Priedmann-Robertson- Walker spacetime, assuming a non-varying a. This will give 
us the usual background behaviour for the usual matter components in the universe: 
pressureless matter, radiation and cosmological constant. In section 1.7 we will then 
introduce the varying-ct model we will use throughout this thesis. 

1.6 The Friedmann- Robertson- Walker Universe 

The standard Friedmann-Robertson- Walker (FRW) cosmological model is based on 
three main theoretical assumptions: The first is that General Relativity is the cor- 
rect description of gravitational interactions, which implies that the model is four 
dimensional. The second assumption concerns the particle content of the universe, 
which is assumed to be described by the Standard Model of Particles (SM). Finally, 
it is based on the Cosmological Principle, which tells us that our universe is spatially 
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maximally symmetric at any constant time and so, isotropic and homogeneous in 
space. These assumptions are based upon three robust observational facts: 

• The observation by Hubble and Slipher that all the galaxies are separating 
from each other, at a rate that is roughly proportional to their separation, is 
realised for a Hubble parameter, H (see later) being approximately constant 
at present, H — Hq > 0. This has been highly verified during the past few 
decades. 

• The relative abundance of the elements with approximately 75% of Hydrogen, 
24% Hehum and other light elements such as Deuterium and Helium-4 with 
small fractions of a percent, is a big success of nucleosynthesis, and at present, 
it is the farthest away in the past that we have been able to compare theory 
and observation. 

• The discovery of the Cosmic Microwave Background Radiation (CMBR) by 
Penzias and Wilson in 1964, signalling the time of photon last scattering. 
This is one of the strongest evidences that our universe started in a stage of 
a hot big bang. The temperature of the CMBR across the sky is remarkably 
uniform: the deviations from isotropy, i.e. differences in the temperature of 
the blackbody spectrum measured in different directions of the sky, are of 
order of 20fiK on large scales, or one part in 10^ [116, 117, 118, 119, 120, 
121, 122, 123, 124, 125, 126]. The observed high degree of isotropy not only 
provides strong evidence for the present level of large-scale isotropy and ho- 
mogeneity of our Hubble volume, but also provides an important probe of 
conditions in the universe at red shifts^ of order 1100. The primeval density 
inhomogeneities necessary to initiate structure formation result in predictable 



^Thc red shift of an object, z, is defined in terms of the ratio of the detected wavelength to the 
emitted wavelength as 

_ Ao a{to) /I fi 1^ 

l + z= — = ——, (1.6.1) 

Ae a[te) 

where a{t) is the scale factor (sec later). Any increase (decrease) in a(t) leads to a red shift (blue 
shift) of the light from distant sources. Since today observed distant galaxies have red shift spectra, 
we can conclude that the universe is expanding. 
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temperature fluctuations in the CMBR, so it can be used to probe theories of 
structure formation. 



Now, a mathematical consequence of the cosmological principle, is that the met- 
ric of the cosmological spacetime, takes the Friedman Robertson Walker (FRW) 
form 

ds"^ = g^^dzf^dx" = df - a^{t)dx^ , (1.6.2) 

where = 0, 1, 2, 3 and dx^ is the metric of the three dimensional maximally sym- 
metric space with constant curvature, ^"^^ R = with parameter k — ±1,0, 
corresponding to a(n) open, closed or flat three slices. This can be written as: 

dr^ 

dx^ = -fijdx'dx^ = — + r^lde^ + sin^ Odcj)^] = dx" + f(x)[dO^ + sin^^ 

1 — kr^ 



where 



and 



(1.6.3) 
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Spatially open, flat and closed universes have different geometries: Light geodesies 
on these universes behave differently, and thus could in principle be distinguished 
observationally. 

The scale factor a{t), can be determined by solving the equations of motion 
coming from an action which contains gravity, and some matter content that we 
describe as a perfect fluid, in consistency with the assumptions of homogeneity and 
isotropy (an observer comoving with the fluid would see the universe around her/him 
as isotropic). The action for this system takes the form: 

J d^xy/^[R + 2A-2K^C^], (1.6.5) 
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where = SttG = Mp^ with G Newton's constant and Mp = 2.4 x lO^^GeV is the 
Planck mass. We have introduced a possible cosmological constant A, and Cm is the 
Lagrangian describing the matter content in the universe. 

The equations of motion derived from (1.6.5) give rise to Einstein's equations 
(where we will assume throughout the thesis h — c= 1), 

G^u = Rnv - ^RQui^ = -i^^T^u - ^Qnv ■ (1.6.6) 
Here the stress-energy momentum tensor for a perfect fluid is defined by 

\/-g 5g^^ 

and can be written as 

= diag(p, -p, -p, -p) , (1.6.8) 

where is the comoving four- velocity, which satisfies, -u^-u^ = — 1 and p{t) and p{t) 
are the energy density and pressure of the fluid, respectively, at a given time in the 
expansion. They are related by the equation of state^: 

p — wp. (1.6.9) 

From Einstein's equations (1.6.6), one obtains the following two (not independent) 
equations of motion. The first one is the Friedmann equation 



H^ = -p-- + -, (1.6.10) 

Qj o 



where if = ^ is the Hubble parameter and a dot means derivative with respect to 
the cosmic time t. the second is the Raychaudhuri equation 

5.4(, + 3rt4. (1.6.11) 



^In general, the parameter w, which gives the speed of sound, can depend on time, but through- 
out this work, we take it to be constant. 
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Now, the ji — Q component of the conservation equation for the stress-energy mo- 
mentum tensor, T^^.^^ — gives the conservation of energy equation: 

p + 3//(p + p) = 0, (1.6.12) 

which is imphed by equations. (1.6.10, 1.6.11). Therefore, after using the equation 
of state (1.6.9) we are left with only two equations that wc can take as the Friedmann 
and the energy conservation, for p(t) and a(t), which can be easily solved. Equation 
(1.6.12) gives immediately 

p ~ a-'(^+'") ; (1.6.13) 
introducing this into Friedmann's equation gives us a solution for a{t) as 

a{t) oc t^(^ (for A = A; = 0) (1.6.14) 

The behaviour of p and a{t) for the typical values of the equation of state are shown 
in Table 1.1 for a flat universe, that is /c = 0. The non flat cases can be found 
straightforwardly. 



Stress Energy 


w 


Energy Density 


Scale Factor a{t) 


Matter 


w = 


p ~ 


a{t) ~ ^2/3 


Radiation 


w = \ 


—4 

p ~ a 


a{t) ~ 


Vacuum (A) 


w = —1 




a{t) ~ exp('sj^t) 



Table 1.1: The behaviour of the scale factor and the energy density for matter, 
radiation and vacuum dominated universes for A; = 0. 



For the more general cases, we can say several things without solving the equa- 
tions explicitly. First of all, note that Raychaudhuri equation (1.6.11) for a vanishing 
cosmological constant A = 0, implies a/a < 0, if {p+3p) is positive. Moreover, a > 
by definition and a/a > since we observe red shifts, not blue shifts. So one can 
conclude that a is a growing function of time at present and that the curve a{t) vs. 
t hits the axis at some finite time in the past, if Einstein's equations are valid at all 
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times. Then i = is defined by a(0) = 0, where the model has a singularity. The 
age of the universe fo, is the time passed since then and has a finite value. 

On the other hand, conservation of energy (1.6.12) requires that p decrease at 
least as fast as as a increases, if p remains positive. Then equation (1.6.10) 
implies that — > a"^ — k (at least) as a ^ oo. Thus the behaviour of the expansion 
depends on the value of A;: for A; = —1, remains positive definite, so a{t) — > i as 
t ^ oo and the universe expands for ever. If A; = 0, d^ also remains positive definite, 
and a{t) f^, with c < 1, as t — > oo. Then the universe expands indefinitely, but 
slower than in the previous case. For A; = 1, the expansion stops at some point, 
d = 0, when 

' - " (1.6.15) 



4 = 



and d takes negative values since d < 0. Then the universe begin to decrease again, 
reaching the singularity in a finite time in the future, see figure 1.1 . 



a(t) A 




Figure 1.1: Three Universes: Closed(K = 1), Flat(K = Q) and Open (K = -1) 

Moreover, using (1.6.10) and (1.6.12) we can write the derivative of the Hubble 
parameter as 

ij = -^(p + p) + A. (1.6.16) 

This equation is already telling us something very important. For flat models. A; = 0, 
this relation implies that reversal from contraction [H < 0) to expansion [H > 0) 
of the scale factor a{t) is impossible since (p + p) > by the weak energy condition 
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(see [127] for a review on the energy conditions in general relativity). 

Let us now introduce a useful concept that will allow us to make general com- 
ments about the cosmological evolution of the FRW universe. This is, the critical 
density, which is defined by: 

Pc^^-1.7x (1.6.17) 

where Hq is the present value of the Hubble parameter, Hq ~ 65Km/s/Mpc~^ ^. 
Using this relation, we can also define the dimensionless density parameter as 

Qi = — , (1.6.18) 

pc 

where the subindex i stands for matter, radiation, cosmological constant and curva- 
ture, today, that is 



A „ 



Using these relations, we can rewrite Friedmann's equation in the form 

l^fc = fi^ + l^A - 1 = VLt,t - 1 (1.6.20) 

where we have neglected Jl^ today. Then, we can write down the following equiva- 
lences: 



k = 


1 




^tot 


> 1 




P> pc 




closed 









^tot 


= 1 




P = Pc 




fiat 




-1 




^tot 


< 1 


44> 


P< pc 




open 



At this point, we have all the information needed to trace the evolution of the 
universe with the assumption that the universe corresponds to an expanding gas of 



*1 Mpc = 3 X 10^^ meters 
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particles described by the Standard Model of particle physics. This gas is assumed to 
be in equilibrium. There are two ways in which the particles leave equilibrium. One 
is when the mass threshold of the particle is reached by the effective temperature of 
the universe, and so, it is easier for the particle to annihilate with its antiparticle, 
rather than being produced again since, as the universe cools down, there is not 
enough energy to produce such a heavy object. The other way is when the interaction 
rate of the relevant reactions F, is smaller than the expansion rate of the universe, 
measured by if. At this time the particles get out of equilibrium. For instance 
at temperatures above 1 MeV the reactions that keep neutrinos in equilibrium are 
faster than the expansion rate but at this temperature H < F and they decouple 
from the hot plasma, leaving then an observable, in principle, trace of the very early 
universe. However, at the present time, we are far from being able to detect such 
radiation (for a detailed analysis, see [52, 128]). 

At a temperature of about lO^K, corresponding to 10 eV, the non-relativistic 
matter content in the universe reached the same density as the rclativistic one, and 
the universe changed from being radiation dominated to being matter dominated. 
When this happened, the expansion rate increased from a(t) oc t^/^ to a(t) oc t^/^. 
From that point on, the temperature decreased more quickly than the (fourth root 
of the) energy density. The first significant event of this epoch occurred at around 
3000K, when decoupling of the radiation from matter occurred. The universe was 
cold enough for atoms to be formed, and so radiation decoupled completely from 
the matter, and the photons propagated freely, giving rise to the cosmic microwave 
background radiation. After this, structure must have formed in the universe, in- 
cluding the first large gravitationally bound systems, such as clusters and galaxies, 
probably due to quantum fluctuations of the early universe, leading to our present 
time. 

1.7 The Bekenstein-Sandvik-Barrow-Magueijo The- 
ory 

In the previous section we have described the evolution of a Friedmann-Robertson- 
Walker universe and its components. Now, in order to study the time and spatial 
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variations of the fine structure constant during the history of the evolution of the 
universe, we need to consider a theoretical model which incorporates variations of 
a into the previous cosmological model. 

Due to the amazing theoretical and observational success of the standard cos- 
mological model, we will choose a varying-a theory which will preserve as much as 
possible the behaviour of an FRW universe. 

As we saw in the previous sections, a common feature, to all theoretical varying- 
a models, is the existence of a scalar field, responsible for the variations of a, which 
is, at least, coupled to the matter fields which interact electromagnetically. Since we 
are interested to investigate the effects of matter inhomogeneities on the evolution of 
the fine structure constant, and we want the preserve Einstein gravity formulation, 
we will then consider a four-dimensional model where the scalar field is only coupled 
to the matter fields which interact electromagnetically. 

With all these desired features in mind, we will choose the Bekcnstcin-Sandvik- 
Barrow-Magueijo (BSBM) theory for a varying fine structure constant. 

Bekenstein did not take into account the effect of the field e in the Einstein 
equations and studied only the time variation of e in a matter dominated universe. 
Sandvik, Barrow and Magueijo have generalised the scalar theory by Bekenstein in 
order to include the gravitational effects of the field, e, responsible for the variations 
of the fine structure constant [107] . In that sense, they have extended the analysis 
by Bekenstein by solving the coupled system of the Einstein equations and the 
Klein-Gordon equation of e [107]. 

In order to simplify Bekenstein's theory, they introduced an auxiliary gauge field 
= eAfj, [107, 88] , and the electromagnetic field tensor can be written now as 

The covariant derivative then becomes — -\- ieoan. To simplify further an- 
other transformation is possible [88]: e — > V = l^e. The field will then be the 
responsible for the variations of the fine structure constant. The scalar field '4> plays 
a similar role to the dilaton in the low energy limit of string and M-theory, with 
the important difference that it couples only to electromagnetic energy. Since the 
dilaton field couples to all the matter (although generally to different sectors with 
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different powers) tfien tfie strong and electroweak cfiarges, as well as particle masses, 
can also vary with the time-position coordinate Xf^. These similarities highlight the 
deep connections between effective fundamental theories in higher dimensions and 
varying-constant theories [58, 79, 78, 69]. 

The total action describing the dynamics of the Universe with a varying-a and 
including the Einstein-Hilbert action for gravity and normal matter takes the form: 

S — J d X\J g {^ilgrav ~l~ ^matter ~l~ ^iji ~l~ ^em^ ) (-^■'^■2) 

where = '^d^tpd^tp, c<j = is a coupling constant, and £em = —\fnuf'^^- The 
gravitational Lagrangian is the usual Cgrav ~ with R the curvature scalar, 

and Chatter is the matter fields Lagrangian. 

To obtain the cosmological equations we vary the action (1.7.2) with respect to 
the metric to give the generalised Einstein equations 

G,. = 87rG (T-/ + + T;-) (1.7.3) 

which are similar to the one derived in equation (1.6.6) but now wc have included 
the energy momentum of ip, which can be obtained using equation (1.6.7). The 
equation of motion for 'ip comes, varying the action (1.7.2) with respect to it 

Uip = -e-2^£em (1-7.4) 

OJ 

The right-hand-side (RHS) of equation (1.7.4) represents a source term for 
which includes all the matter fields which are coupled to it. These include not only 
relativistic matter (like photons), but as well as non-relativistic one that interact 
electromagnetically. 

It is clear that Cem, vanishes for a sea of pure radiation since £em = (-^^ — B"^) /2 = 
0. The only significant contribution to a variation of ip comes from nearly pure elec- 
trostatic or magnctostatic energy associated to non-relativistic particles. In order to 
make quantitative predictions we then need to know how much of the non-relativistic 
matter contributes to the right-hand-side (RHS) of equation (1.7.4). This can be 
parametrised by the ratio C, = Cgm/ Pm, where Pm is the energy density of the non- 
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relativistic matter [129]. 

For protons and neutrons, ( can be estimated from the electromagnetic cor- 
rections to the nucleon mass, 0.63MeV and —O.lSMeV, respectively [56]. This 
correction contains the contribution, which is always positive, and also terms 
of the form j^a^, where is the quarks' current [108]. Hence we take a guidance 
value of C ~ 10~^ for protons and neutrons. 

Using the parameter ^, the fraction of electric and magnetic energies may then 
be written as: 

= - C = - (1-7.5) 

Pm Pm 

where E'^ and B'^ are the electric and magnetic energies respectively. Using equation 
(1.7.5), then equation (1.7.4) becomes 

Di^ = -e-'^Pm (C^ - C^) (1-7.6) 

CO 

Since we are interested in the cosmological evolution of a, instead of using both 
parameters and C^, we will use throughout this thesis, the cosmological param- 
eter, C, defined as C = ~ C^j which in the limit where ^ is positive, and 
when <^ is negative. 

Note that, the cosmological value of ( has to be weighted, not only by the 
electromagnetic-interacting baryon fraction, but also by the fraction of matter that 
is non-baryonic, for instance dark matter. The coupling to dark matter is motivated 
by two aspects. The first is the fact that dark matter might be electromagnetically 
charged, for instance superconducting cosmic strings, and if that is the case the scaler 
field is necessarily coupled to it. The second is the Olive and Pospelov generalisation 
of the Bekenstein model, where they claim that the field responsible for the variations 
in a may be more strongly coupled to dark matter than to baryons (depending on 
the nature of dark matter) [57]. Hence the value and sign of ( depends also on the 
nature of dark matter to which the fields ip might be coupled. For instance, ( = —1 
for superconducting cosmic strings, where Cem ~ —B'^/2 [108, 109], and in the case 
of neutrinos |^| <^ 1. 



The universe, we will be studying, will be described by a flat, homogeneous and 
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isotropic Priedmann metric. The universe contains pressure-free matter, of density 
Pm-i a cosmological constant A, of density px = A/{8nG) and radiation, of density 
Pr- The Priedmann equation can be obtained from the Einstein equations (1.7.3), 

H'^^ {pm (1 + Id e-'^) + pre-'^ + p^ + pa) (1.7.7) 

where = ^ip"^. Using the cosmological parameter (, the evolution equation for -0 
comes from 1.7.4 

xi^ + 3HiP^ -- e-'^CPm- (1.7.8) 

The conservation equations for the non-interacting radiation and matter densities, 
Pr and pm respectively, are: 

p'm + SHpm = 
Pr + 4:Hpr = 2V'Pr (1.7.9) 

so Pm oc a~^. The last relation in equation (1.7.9) can be written as 

%+AHpr^O, (1.7.10) 

with Pr = p^e"^^ oc a~^. 

Equation (1.7.8) may be expressed in terms of the kinetic energy density of the 
ip field, p^ = to give 

p'^ + 6Hp^ = 2\ -e'^'^CmPm^/P^- (1.7.11) 
V CO 

The ip field behaves like a stiff Zeldovich fiuid with p^ oc when the RHS vanishes. 

The Bekenstein-Sandvik-Barrow-Magueijo theory, enables the cosmological con- 
sequences of a varying fine structure constant to be analysed self-consistently. Equa- 
tions (1.7.7), (1.7.8 and (1.7.9), govern the Friedman universe with time-varying a . 
These equation can be evolved numerically from early radiation epoch, through the 
matter era and into vacuum domination. Notice from Figure 1.2 and 1.3 that the 
energy density associated to the scalar field ip is always negligible at all the stages 
of the universe history. 
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Figure 1.3: Evolution of log Pm (solid lines) , logp^ (dash-dotted line), logp;^ (dotted 
line) and log pr (dashed line) as a function oflog{a). 
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The inclusion of a varying-o; will then affect very little the cosmological expansion 
of the universe. Hence, the Hubble diagram will be precisely the same as that of a 
universe, with a constant a, with ^ 0.3, ^ 0.7 and h — 0.65 [107]. This is 
the effect that we desired in a varying-a model: the less possible deviations from 
the standard FRW universe. 



CHAPTER 2 



Background Solutions for the Fine Structure 

Constant 



Finjo viver num tempo sem fim 
e de tanto fingir. . . 
esqueci-me que so sei viver assim 
- Ernesto Mota - 



2.1 Introduction 

In the last chapter we have described the BSBM varying-o; cosmological model we 
will use to study the effect of the matter inhomogeneities in the cosmological evo- 
lution of the fine structure constant. Before starting to investigate those effects, we 
first need to study the evolution of a in the homogeneous and isotropic background 
universe and to calculate analytical solutions of the equation of motion (1.7.8) for ■0. 
These analytical solutions will be useful when we study the existence and behaviour 
of the inhomogeneities. 

We will then provide a complete analysis of the behaviour of the solutions of the 
non- linear propagation equation (1.7.8) for appropriate behaviours of a{t), for each 
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epoch the universe went through: radiation dominated, dust dominated and dark 
energy dominated. 

The aim of this chapter brings a question: How shall we solve the non-linear 
differential equation (1.7.8)? One cannot hope to obtain exact solutions to most 
non-linear differential equations. The reason is that there are only a limited number 
of systematic procedures for solving them, and these apply only to a very restricted 
class of equations. Moreover, even when a closed-form solution is known, it may be 
so complicated that its qualitative properties are obscure. Thus, for most non-linear 
equations it is necessary to have reliable techniques to determine the approximate 
behaviour of the solutions. In order to calculate and to understand the exact solu- 
tions of equation (1.7.8), we will use a phase-space analysis. Prom this analysis, we 
will be able to find stable asymptotic solution which can be considered to describe 
the behaviour of -0 during the different eras the universe went through. 

2.2 An Approximation Method 

We have seen in the last chapter, that the energy density of the field associated to 
variations of a is always negligible with respect to all the other energy content of 
the universe (p^ <^ Pm, Pr, Pa)- This means the Friedmann models with varying a 
have the property that the homogeneous motion of the ip does not create significant 
metric perturbations at late times [108]. So, far from the initial singularity, we can 
safely assume that the expansion scale factor follows the form for the Friedmann 
universe containing the same fiuid when a does not vary [( = = ip). The behaviour 
of ijj then follows from a solution of equation (1.7.8) in which a{t) has the form for a 
Friedmann universe for matter with the same equation of state in general relativity 
when ( — — t/j. This behaviour is natural. We would not expect that very small 
variations of the coupling to electromagnetically interacting matter would have large 
gravitational effects upon the expansion of the universe. Thus, in this chapter we 
will provide a complete analysis of the behaviour of the solutions of equation (1.7.8) 
for appropriate behaviours of a{t). 

Before starting the analysis and to calculate the asymptotic exact solutions let 
us show that the assumptions we claim are indeed right. Hence, we will show how 
good is the approximation of assuming that the scale factor, a{t), will not deviate 
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from the the usual FRW universe that we have calculated in the previous chapter, 
see table 1.1. 

We consider spatially flat universes {k = 0) and assume that the expansion scale 
factor is that of the Friedmann model containing a perfect fluid: 

a = r (2.2.1) 

where n is a constant. The late stages of an open universe containing fluid with 
density p and pressure p obeying p + 3p > can be studied by considering the case 
n—1. We rewrite the wave equation (1.7.8) as 

(-00^) = p^a^ exp[-2'0] 

(jj 

Therefore, since PmO^ is constant this reduces to a Liouville equation of the form 

(ipa^y^ A^exp[-2V'] (2.2.2) 
where N is a constant, defined by 

N = -'^Pma^ ■ (2.2.3) 

We shall consider first the cosmological models that arise when the defining 
constant N is negative. This arises when the constant C < 0; indicating that the 
matter content of the universe is dominated by magnetic rather than electrostatic 
energy. The value of ( for baryonic and dark matter has been disputed [56, 57, 107]. 
It is the difference between the percentage of mass in electrostatic and magnetostatic 
forms. As explained in the previous chapter and [107] , we can at most estimate this 
quantity for neutrons and protons, with Cn ~ Cp ~ 10~^. We may expect that 
for baryonic matter (. ~ 10~^, with composition-dependent variations of the same 
order. The value of C for the dark matter, for all we know, could be anything 
between -1 and 1. Superconducting cosmic strings, or magnetic monopoles, display 
a negative unlike more conventional dark matter. It is clear that the only way to 
obtain a cosmologically increasing a in BSBM is with C < 0, i.e with unusual dark 
matter, in which magnetic energy dominates over electrostatic energy. In [107] it 
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was showed that fitting the Webb et al results requires C,/u = — 2 ± 1 x 10~^, where 
C is weighted by the necessary fractions of dark and baryonic matter required by 
observations of the gravitational effects of dark matter and the calculations of Big 
Bang nucleosynthesis. We note also that in practise C, might display a significant 
spatial variation because of the change in the nature of the dominant form of matter 
over different length scales. For example, a magnetically dominated form of dark 
matter might contribute a negative value of C, on large scales while domination of 
the matter content by baryons on small scales would lead to C > locally. We will 
not investigate the effects of such variations in this thesis. 

2.2.1 The validity of the approximation 

We have assumed that the scale factor is given by the FRW model and then solved 
the evolution equation. This is a good approximation up to logarithmic correc- 
tions. Here is what happens to higher order. 

We take the leading order behaviour in (1.7.7) 

(^)'«^ YP^(l + IC|e-^^) (2.2.4) 

Now if we take a = t^^^ so 

Stt _ 4 

and solve equation (2.2.2) we get asymptotically, 

= - \n[2Nhit] ln[lni] (2.2.5) 

to leading order at late times. Suppose we now re-solve (2.2.4) with the |C|e~^^ 
correction included 

Note that the kinetic term which we neglected is of order 
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and so is smaller than the term we have retained. Solving equation (2.2.4) we have 

2|Cl 



Ina = 



'Int 



Int 



Int 



sh-H )■ + W^W1 + ^ 



ICI 



ICI 



ICI 



(2.2.8) 



Note that when |C| — >• this gives the usual Ina = |lnt. When |d 7^ we have 



In a 



2 ICI 



and 



a = t2/3(int) 



ICI/3 



(2.2.9) 



where |C| is small and so the corrections to the a — t^l^ ansatz are small. In terms 
of the Hubble rate: 

H _ 2 ICI 
3i 3Uni 

If we include the kinetic corrections to equation (2.2.4) then as -0 — ^ 

2 A 



2tlnt 



^2 1 + 



9^2 



+ 



S 



ICI 



(2.2.10) 



where 



So, ii X = Int, we have 

y/x-'+\C\x + S + S J 



3, 

-Ina 



dx 



XyJ x^ + ICI X + 5" 



dx 



V^2 + ICk + ^ 



(2.2.11) 



3, 

-Ina 



= Vx2 + |C|x + >5-ln 



25 + Id X + 2V^V^2^W+^' 



2^x2 + |C|a; + >S + 2x+ Id 



(2.2.12) 



Again, as i — > 00 the leading order behaviour is that found in equation (2.2.9). 
In the radiation era, where a{\) oc f^/^, we have an exact solution of equation 
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(2.2.13) 



so the corrections to the Friedmann equation look hke 




2 



1 

4i2 




(l + |C|e-^) 



(2.2.14) 



and these corrections fall off much faster than in the dust case. Again, our basic 
approximation method holds good to high accuracy. 

2.2.2 A linearisation instability 

Despite the robustness of the basic test-motion approximation that we are employing 
to analyse the evolution of ilj{t) as the universe expands, there is a subtle feature 
the non-linear evolution equation (2.2.2) which must be noted in order that spurious 
conclusions arc not drawn from an approximate analysis. We sec that the right- 
hand side of equation (2.2.2 ) is always positive. Therefore ip can never experience 
an expansion maximum (where ip = and ip < 0) and therefore ip{t) can never 
oscillate. However, if we were to linearise equation (2.2.2), obtaining 



then for > 1/2 the right-hand side takes negative values and pseudo-oscillatory 
solutions for ip would appear that are not the linearised approximation to any true 
solution of the non- linear equation ( 2.2.2). Care must therefore be taken to ensure 
that analytic approximations are not extended to large ip and that numerical anal- 
yses arc not creating spurious spirals in the phase plane by virtue of a linearisation 
procedure; for a fuller discussion sec rcf. [110]. 

These considerations can be taken further. It is possible for ilj{t) to decrease, 
reach a minimum and then increase. But it is not possible for ilj{t) to decrease if 
it has ever increased. A second interesting consequence of this feature of equation 
(2.2.2) is that it holds true even if a{t) reaches an expansion maximum and begins to 
contract. Thus in a closed universe we expect ip and a to continue to increase slowly 
even after the universe begins to contract. This will have important consequences 



(ipa^y ^ A^exp[-2V'] A^(1-2V' + 0(V'^)) 
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for the expected variation of ijj and a in realistically inhomogeneous universes. 



2.3 Two-Dimensional Non-Linear Autonomous Sys- 
tems 



Autonomous systems of equations, when they are interpreted as describing the mo- 
tion of a point in the phase space, are particularly susceptible to some very beautiful 
techniques of local analysis. By performing a local analysis of the system near what 
are know as critical points, one can make remarkably accurate predictions about the 
global behaviour of the solution. 

In this section, we shall not attempted to give a proper description of the the 
dynamical systems field, but we will only state the basic results needed to understand 
this chapter. For a better and proper study of the topic see for instance [130, 131, 



Differential equations which do not contain the independent variable explicitly 
are said to be autonomous. Any differential equation is equivalent to an autonomous 
equation of one higher order. 

It is convenient to study the approximate behaviour of an autonomous equation 
of order 2 when it is in the form of a system of 2 coupled first-order differential 
equations. The general form of such a system is 



where the dots indicate a differentiation with respect to the independent variable, 
for instance t. 

The solution of the system (2.3.1) is a curved or trajectory in a two-dimensional 
space called phase space. The trajectory is parametrised in terms of t : yi — yi{t), 
y2 = 1/2 (t). 

We assume that /i, /2 are continuous differentiable with respect to each of their 
arguments. Thus, by the existence and uniqueness theorem of differential equations 
[130] any initial condition yi{0) = ai, ^2(0) = 02, gives rise to a unique trajectory 



132]. 



yi = /i(yi,y2) 
y2 = /2 (1/1, 1/2) 
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through the point (01,02). To understand this uniqueness property geometrically, 
note that every point on the trajectory [yi{t),y2{t)], the system (2.3.1) assigns a 
unique velocity vector [yi{t),y2{t)\ which is tangent to the trajectory at that point. 
It immediately follows that two trajectories cannot cross; otherwise, the tangent 
vector at the crossing point would not be unique. 

2.3.1 Critical points in phase space 

If there are any solutions to the set of simultaneous algebraic equations 



then there are special degenerate trajectories in phase space which are just points. 
The velocity at these points is zero so the position vector does not move. These 
points are called critical points. 

Studying the phase space it is possible to make elegant global analyses of the 
system. The possible behaviours of a trajectory in a two dimensional system [130] 
are: 

• The trajectory may approach a critical point as i — > +00. 

• The trajectory may approach 00 as i ^ +00. 

• The trajectory may remain motionless at a critical point for all t. 

• The trajectory may describe a closed orbit or a cycle. 

• The trajectory may approach a closed orbit (by spiralling inward or outward 
toward the orbit) as i — > +00. 

The possible local behaviours for trajectories near a critical point are: 

• 1. All trajectories may approach the critical point along curves which are 
asymptotically straight lines as i — > +00. We call such a critical point a stable 



/i(yi,y2) = 

/2(yi,y2) = 



(^ yi = 0) 
(^ m = 0) 



(2.3.2) 



node. 
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• 2. All trajectories may approach the critical point along spiral curves as i — > 
+00. Such a critical point is called a stable spiral point. 

• 3. All time reversed trajectories, that is, y{t) with t decreasing, may move 
toward the critical point along paths which are asymptotically straight lines 
as t ^ —oo. Such a critical point is an unstable node. As t increases, all 
trajectories that start near an unstable node move away from the node along 
paths that arc approximate straight lines, at least until the trajectory gets far 
from the node. 

• 4. All time-reversed trajectories may move forward the critical point along 
spiral curves as i — > — oo. Such a critical point is called an unstable spiral 
point. As t increases, all trajectories move away from an unstable spiral point 
along trajectories that are, at least initially, spiral shaped. 

• 5. Some trajectories may approach the critical point while others move away 
from it as t — > +00. Such a critical point is called a saddle point. 

• 6. All trajectories may form a closed orbit about the critical point. Such a 
critical point is called a centre. 

2.3.2 Matrix methods 
Linear autonomous systems 

Since two-dimensional linear autonomous systems can exhibit any of the critical 
point behaviours that we have described above, it is appropriate to study linear 
systems before going to non-linear ones. With this in mind we introduce an easy 
method for solving linear autonomous systems. 

A two dimensional linear autonomous system yi = ayi + by2, y2 = cyi + dy2 may 
be re- written in matrix form as 

Y = MY (2.3.3) 



where M — 




2.3 Two-Dimensional Non-Linear Autonomous Systems 



38 



It is easy to verify that if the eigenvalues Ai and A2 of the matrix M are distinct 
and Vi and V2 are eigenvectors of M associated to the eigenvalues Ai and A2, then 
the general solution to equation (2.3.3) has the form 

Y{t) = CiVie^i* + CsVae^^* (2.3.4) 

where Ci and C2 are constants of integration which are determined by the initial 
position Y(0). 

The hnear system (2.3.3) has a critical point at the origin (0,0). It is easy to 
classify this critical point once Ai and A2 are known. Note that, Ai and A2 satisfy 
the eigenvalue condition 

det[M - IX]^det \"'~^ ^ | = A^ - A(a + d) + ad - 6c = 
y c d — X I 

If Ai and A2 are real and negative, then all trajectories approach the origin as 
f — > 00 and (0, 0) is a stable node. Conversely, if Ai and A2 are real and positive, 
then all trajectories move away from (0, 0) as i — > 00 and (0, 0) is an unstable node. 
Also, if Ai and A2 are real but Ai is positive and A2 is negative, then (0, 0) is a saddle 
point; trajectories approach the origin in the V2 direction and move away from the 
origin in the direction Vi. 

Solutions Ai and A2 of (2.3.2) may be complex. However, when the matrix M 
is real, then Ai and A2 must be a complex conjugate pair. If Ai and A2 are pure 
imaginary, then the vector Y(t) represents a closed orbit for any Ci and C2 and the 
critical point at (0, 0) is a centre. If Ai and A2 are complex with non-zero real part, 
then the critical point at (0, 0) is a spiral. When Re Ai,2 < 0, then Y(t) 0, (0, 0) 
is a stable spiral point; conversely, when Re Ai_2 > 0, (0, 0) is an unstable spiral 
point. 

Two-dimensional non-linecir systems 

The analysis of a non-linear system is a much more complicate case. Nevertheless 
non-linear system can be analysed near its critical points, and the result of that 
study can be indeed very accurate if the system is almost linear. In our particular 
case, the condition for almost-linearity can be checked comparing our results with 
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the numerical integrations performed in [108]. We will see that our results are indeed 
very good approximation to the solutions found numerically, with exception to one 
case where one of the eigenvalue solution of equation (2.3.2) is zero. The stabihty 
analysis of this case will be described bellow in 2.3.2. 

The local analysis of a non-linear system, to which the linear approximation 
works well, consist into linearise the equations and then procedure as usually done 
for a linear system. We first identify the critical points. Then we perform, a local 
analysis of the system very near them. Using matrix methods, we identify the nature 
of the critical points of the linear system. Finally, we assemble the results of our local 
analysis and synthesise a qualitative global picture of the solution to the non-linear 
system. 

Critical point with one zero eigenvalue 

In this section we will only present the results obtained in [133] and we refer the 
reader to that reference for a detailed explanation. 

If one of the eigenvalue solutions, Ai,2, of equation (2.3.2) has, for instance, 
Re{Xi) — and Re{X2) < then the stabihty cannot be decided by the linear 
approximation. 

In general, a non-linear autonomous system can be written in the form 

Y = MY + 7V(Y) (2.3.5) 

where M is a constant matrix 2x2 and N{Y) is non-linear in Y. The linearisation 
of this system gives (2.3.3). Consider now that, the two eigenvalues of M, have 
it!e(Ai) = and Re{X2) < 0, corresponding to the eigenvectors Vi and V2. In order 
to perform the stability analysis of the critical point (0, 0), we should procedure as 
follows (For other alternatives to this method see [131]). 

• First, we apply a linear transformation to Y using the matrix formed by the 
eigenvectors Vi and V2. This will split the system into critical and non-critical 
variables; the critical variable is the eigenvector corresponding to i?e(Ai) = 0, 
in our case Vi. Hence, Y — > X = {xi,X2). In the new coordinates, we now 
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have that the following system 



xi = ni{xi,X2) 

X2 = P2lXl+P2Xl+n2{Xi,X2) 



(2.3.6) 



where ni and n2 are at least quadratic in both Xi and X2- 

• Secondly, if there are linear terms in xi, in the above system, then they must 
be eliminated. That can be achieved solving the equation X2 — 0, which will 
give us X2 has a function of xi (the critical variable). Say, for instance, that 
solution give X2{xi). 

• At last, we perform the following transformation {xi, X2) {zi, Z2) defined as 



Having performed this steps we end up with an autonomous system for the 
coordinates (-21,-22) 



where the functions hi and /i2 are at least quadratic in their arguments, the stability 
of the critical point Y = (0, 0) is determined by the stability of Z = {zi, Z2) — (0, 0). 

Since the non-null eigenvalue of M has a negative real part {Re{\2) < 0, by 
assumption), the coordinate Z2 is asymptotically stable about the origin (0,0) as 
t — > 00. The asymptotic stabihty of the critical variable is therefore determined as 
t — > 00 and -^2 — by the leading term of hi{zi, 0). It can be shown [133] that when 
the first power of hi{zi,Q) is even (recall that hi is at least a quadratic function), 
then (-21,-22) = (0,0) is unstable. 



{zi,Z2) = {xi,X2 - X*2 {xi)) 



(2.3.7) 



ii = hi{zi,Z2) 

Z2 = P21X1 + h2{Zi, Z2) 



(2.3.8) 
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2.4 Phase-Plane Analysis 

We are ready now to analyse the phase space of equation (2.2.2). We first transform 
it into an autonomous system of two first order differential equations. Then we 
identify the critical points and we perform, a local analysis of the system very near 
them. The exact system will be well approximated by a linear autonomous system 
near the critical points. This will be seen when we compare our results with the 
numerical results performed in [108]. Using matrix methods we identify the nature 
of the critical points of the linear system. Finally, we assemble the results of our local 
analysis and synthesise a qualitative global picture of the solution to the non-linear 
system. 

2.4.1 A transformation of variables 

In this section we will look at the equation of motion (2.2.2 ) when the expansion 
scale factor takes the power-law form (2.2.1). The evolution equation for the field 
then becomes: 

^ (^jjjt^^^ = Nex.p[-2^] (2.4.1) 
with N > 0. We introduce the following variables : 

3n 1 

x^ln{t) y^ip-Ax-B A^l-— B ^ -ln{N) (2.4.2) 



and rewrite (2.2.2) as: 



y" + (3n-l)(y' + l-^) = e-^^ (2.4.3) 



where ' = d/dx. This second-order differential equation can be transformed into an 
autonomous system by defining v — ^ and u — y: 

^ ^ e-^" + (l-3n)(^ + l-^) (2.4.4) 

du 
dx 



With the aim to obtain asymptotic exact solution, we have to search for the 
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critical points of equation (2.4.4). The critical points of a system are defined as 
the values of u and v that give v' — Q and u' — 0, where the prime represents a 
derivative with respect to x. The importance of the critical points is the fact that 
we can linearise the system around then, and if they are stable, the solution of the 
system near them, can be used as an asymptotically stable solution. 

Looking at equation (2.4.4), it is straight forward to see that the system has a 
finite critical point at 

(x.e,^c) = (-2ln[(l-3n)(--l)],0) (2.4.5) 

when n G (|, |) and it has a infinite critical point at (-Uc, v^) = (+oo, 0) when n = ^ 
or n = |. The finite critical points correspond to the family of exact solutions of 
the form ip = B + Aln{t) found in [108]. In the original variables these solutions are 

, 1 , , 2N , /2 -3n\ , , , 

In order to analyse the system fully, we will study finite critical points and the 
critical points at infinity separately. We will also distinguish several domains of 
behaviour for n: n = |, n = |, n G (|, |), n < | and n > ^. In some cases we will 
supplement our investigations by numerical study of the system (2.4.4). 

2.4.2 The finite critical points 

The cases n G (|, |) 

In the domain where n G (|, |) the system (2.4.4) does not appear to be exactly 
integrable. So in this section we will study the behaviour of the system near the 
critical points and at late times. 

There is a finite critical point {uc,Vc) — (— | In [(1 — 3n)(^ — 1)], 0) for n G 
(|, |). Linearising the system 2.4.4 about it we obtain: 

dV 

— = (3n-l)(3n-2)C/+(l-3n)y (2.4.6) 
ax 

dU 
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where V — v — Vc and U — u — Uc- Since the characteristic matrix is non singular 
the critical point is simple. Hence the non-linearised version of this system has the 
same phase portrait at the neighbourhood of the critical point. 

The eigenvalues ^1^2 and corresponding eigenvectors x^^'^ of the system (2.4.6) 
are: 



. 1 - 3n ± V9 - 42n + 45n2 

U,2 - ^ (2.4.7j 

, 1 - 3n ± V9 - 42n + 45^2 , 
X*''' = ( ^ ,1) 

These eigenvalues will be complex numbers for n e (|, |) and they will be pure real 
numbers when n e (|, |). Notice that for both cases the real part of the eigenvalues 
is always negative, so the critical point is a stable attractor. The general solution 
of the linearised system (2.4.6) can be expressed as: 

{U{x), V{x)) = Ae^^^'xi + 5e«^^X2 

where A, B are arbitrary constants and x = ln{t). From the transformations (2.2.2) 
we can obtain explicitly the expression for ilj{t) near the critical point. There are 
three possible behaviours of the solutions near the critical point: 

Pseudo-oscillatory t/j behaviour 

When n G (|, |) the linearised evolution for the ip field exhibits damped oscillations 
about its asymptotic solution as t ^ 00: 



where 



+1" (Bcos[/31n(t)] +^sin[/31n(t)]) 



a = (2.4.9) 
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We note that this family of solutions and asymptotes includes the case of the 
radiation-dominated universe {n — 1/2). The phase portrait and the ip vs. Int 
for this case are shown in figure 2.1. But in the case of a universe containing the 
balance of matter and radiation, displayed by our own, it need not be the case that 
the asymptotic behaviour, displayed by the exact solution for the critical point, is 
reached before the radiation-dominated expansion is replaced by matter-domination, 
see references [107] and [108] for further discussion of this point. However as we dis- 
cussed above these oscillations are an artifact of the linearisation process and the 
part of the solution ( 2.4.8) controlled by the constants A and B is only valid for 
small times, hence we called this behaviour pseudo-oscillatory. 



1 ^ ^ -130 




-0.2 0.2 



Figure 2.1: Numerical plots of the phase space in the {w,v) coordinates, and the ip 
evolution with x — log{t) for n — ^ . The '+ ' sign is a saddle point and the triangle 
is a stable node. 



Non-oscillatory behaviour 

When n G (|,|), the ip field will approach its asymptotic behaviour in a non- 
oscillatory fashion as t — > cxo: 

= i'"' (l-3„H3,.-2) ' P-^-"' 

+1- (2 - 3n) \n{t) + t^'^ {A + B^^) 
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where 



S = 



3n 



7 



42n + 45n2 



(2.4.12) 
(2.4.13) 



Note that 5 + 7 < in the domain n E |) , so at late times {t — > oo), both the 
pseudo-oscillatory and non-oscillatory solutions case will approach the asymptotic 
solution defined by the appropriate value of n. 

Intermediate behaviour 

A transition between the pseudo-oscillatory and non-oscillatory regimes happens 
when P — and JT- = | ; the t/j field then has the following solution in the vicinity of 
the critical point: 



V^(i) = iln(7V) + lln(i)+^rt. 



(2.4.14) 



The phase plane structure and evolution of ijj vs. Int for this case are shown in 
figure 2.2. 




Figure 2.2: Numerical plots of the phase space in the {w,v) coordinates, and the ip 
evolution with x ~ log{t) for n—^. The '+ ' sign is a saddle point and the triangle 
is a stable node. 
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Overview 

The late-time solution of the ip field solution is given by asymptotic behaviour of 
equations (2.4.8) or (2.4.11), which is: 

^-i'"' (l-3n)(3„-2) '4p-^"""W 

This shows that the solutions (2.4.8) and (2.4.11) generalise the ones found in [108], 
since they can be obtained setting A — B — {)'m.{ 2.4.8). In particular, the case of 
a radiation-dominated Universe {n — |), we have from (2.4.8): 

m = ^(log(8)+log(7V)+log(i)) (2.4.15) 

. , „ ,5V31og(t), , . .5V^log(t), 
+^4 ^cos[ ^ ] + Asin[ /'^ 



A full mathematical summary of the change in structure of the phase space 
with changing n for the system (2.2.2) is given in section 2.6. There wc include 
cosmologically unphysical values of n and show how the critical point structure 
bifurcates with the change in value of n. 

2.4.3 The critical points at infinity 

In order to describe the qualitative evolution of the system, we must determine the 
behaviour of the system (2.4.4) near the critical point {uc,v^ — (+00,0). In order 
to bring the critical point to a finite value, we define: 

Using the new coordinate w we can re- write the system (2.4.4) as: 

— = -2wv 2.4.16 
ax 



: (l-3n) f 1-y + -{;) 



dx 



This system has critical points on the {w, v) plane when (0, — 1+^^) and ((1 — 3n) ( 



3n 
2 



2.4 Phase-Plane Analysis 



Note that the second critical point is just the same as the one we have analysed in 
the previous section. In this subsection we will then only analyse the critical point 
{wc, Vc) — (0, —1 + ^) since it corresponds to the case where u — > +oo. 

Proceeding as before, and linearising (2.4.16) about {wc,Vc) — (0,-1 + ^) we 
obtain: 

dW 



dx 
dV 

dx 



(2 - 3n)W (2.4.17) 
W + {1- 3n)V 



where V — v — Vc and W — w — Wc- Again, the characteristic matrix of the system 
is non singular, and the critical point is simple. Hence, system (2.4.17) will have 
the same phase portrait as (2.4.16) in the neighbourhood of the critical point. 

The eigenvalues ^1^2 and corresponding eigenvectors x^^'^ of the system (2.4.17) 
are: 

6 = l-3n, X^^ = (0,1) (2.4.18) 

6 = 2-3n, X^^ = (l,l) (2.4.19) 

These eigenvalues are always real. The critical point is an attractive node forn > |, 
a saddle point when n e (|, |), and it will be an unstable point when n < ^. The 
general solution of the system (2.4.16) in the neighbourhood of the critical point 
{wc, Vc) = (0, — 1 + ^), can be expressed as: 

{w{x) - Wc, v{x) - Vc) = Ae^^'^xi + Be^^^X2 
where A, B are arbitrary constants and x — ln(t). Therefore, near the critical point 

{Uc,Vc) = ( + 00,0): 

m = lH^) (2.4.20) 

so the scalar field is constant. 

Note that in the domain n G (|) |) this is just a transitory solution since the 
critical point is a saddle point. In the domain n < | it is an unstable critical 
point, possibly relevant as an early-time solution to braneworld cosmologies in the 
high-density regime where the Hubble expansion rate of the universe is linearly 
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proportional to the density, so n = 1/4 for a braneworld containing radiation, n — 
1/12 for a massless scalar field, and n = 1/6 for dust. However, in the i — > limit 
the assumption that the ■^^ and Cexp[— 2'0] terms can be neglected in the Priedmann 
equation (1.7.7) will break down. 

The cases of n > 2/3 and de Sitter expansion 

An interesting case exists when tt- > |, since the critical point is a stable attractor 
and so this means that the constant-^' behaviour (2.4.20) is the late-time attractor, 
in agreement with the conclusions of references [107] and [108]. This is an important 
feature of a universe which exhibits accelerated expansion in its late stages {n> 1). 
It means that the present value of a is the asymptotic one. It also means that 
variations of a are turned off by the domination of the expansion dynamics by 
negative curvature or by any quintessence field [107], [108]. This property may 
provide important clues to explaining why our universe possesses small but finite 
curvature or quintessence energy today: if it did not then the fine structure constant 
would continue to increase until it was impossible for atoms and molecules to exist 
[109] 

In the n > 2/3 case we can find a detailed asymptotic solution for equation ( 
2.4.1) which has the form 

= c + 5(t + t,f-'- + ^'^^}~l^\ t + to)^-'" 

with C, to constants. We see immediately that for this solution -0 approaches a 
constant as i — > oo. Of particular interest is the case of a curvature-dominated open 
universe, which has n — 1. The phase plane structure and the evolution of vs. 
\n.t for this case are shown in figure 2.3. 

We see that this approach to constant behaviour occurs for all universes that 
accelerate {n > 1) and so we would expect to find it also in the case of a de Sitter 
background universe with 

a{t) — exp[Ai] 

where A > is constant. 
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Figure 2.3: Numerical plots of the phase space in the {w,v) coordinates, and the ip 
evolution with x — log{t) for n — 1. The '+ ' sign is a saddle point and the triangle 
is a stable node. 

Substituting this in equation (2.2.2) we find a late-time asymptotic solution 
^ = C + Bexpl-SXt] - ^^^^^ ^ cxp[-2C - 3Xt] C 



as i — > oo. 

Notice that we were not able to fully describe the nature of this critical point in 
the cases where n = | or n = | since one of the eigenvalues of the systems is zero. 
In order to do so, we will study these two cases individually, in particular the case 
n = I is important since it describes the scale factor evolution on a dust dominated 
universe. 



The n = ^, a{t) — ts case: an exact solution 

The n = ^ case can be exactly integrated. It is of interest as an exact solution in 
its own right but it corresponds to the case of a universe whose expansion dynamics 
are dominated by the effects of a fluid with a p = p equation of state, or a massless 
scalar field. It also describes the behaviour of the i/j field in an anisotropic universe 
of the simple Kasner type. 
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We see when n — 1/3 the system (2.4.4) has the form: 

^ = e-^" (2.4.21) 
ax 

du 
dx 

This integrates to give 

where is a constant. Hence, we have two possible solution branches: v{u) = 
E — e~^". Both, positive (+) and negative (— ) solutions for v lead to the same 
result. Choosing the positive branch of the v{u) solution and using equations (2.4.21) 
we obtain: 

u(x) = - \n(E) + ln[cosh {E {Ci - x))] 
where Ci is an integration constant. From (2.4.2) we have a solution for ip{t): 

In(^) 1 
jjj{t) = - ln{E) + - ]n{t) + ln[exp[SCi]r^ + t^exp[-^Ci]] 

and the asymptotic limit when i — > oo gives 

^l;{t)^i^ + E)\n(t + to). 
with to constant. The phase space and ip{t) vs. Int evolution is shown in figure 2.4. 



The n = I , a = ts case 

The case of a dust-dominated universe is mathematically special because of the 
presence of a zero eigenvalue in the stability analysis. It is illuminating to consider 
this case separately with an asymptotic analysis that extends the earlier study in 
[108]. 

Consider equation (2.4.1) with n = 2/3. If we introduce the new variable x — 
ln(t), then 

^" + V^' = iV exp[-2^] (2.4.22) 
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Figure 2.4: Numerical plots of the phase space in the {w,v) coordinates and time 
evolution of ip^x) n = |. The central line which passes through the origin is an 
attractor for w > and an unstable line for w < 0. 



At large x, the asymptotic form of this equation has the form: 

= - H2N{x + xo)] - - J] ~ , + C exp[-x] (2.4.23) 

where C and xq> Q are arbitrary constants. The leading order behaviour as t — > oo 
is therefore (cf. (2.2.5)) 

~ -ln[2Arin(i)] 
2 

Stability analysis of the n = 2/3 asymptote : 

Performing the coordinate transformation u = —^ln[w) on the system (2.4.4) 
when n = 2/3, we have: 

dv , ^ ,v 

— ^ w-v (2.4.24) 

dx 

dw 

— —2vw 

dx 

This system has a critical point at the origin of the {w, v) plane. It corresponds to the 
case where m — > oo. The linearisation about this critical point gives two eigenvalues: 



2.5 Some General Asymptotic Features 



52 



Ai = —1 with the eigenvector xi = (0; 1) A2 = with the eigenvector X2 — 
(1,1). Since we have a zero eigenvalue, the stabihty is determined by the non- 
hnear behaviour and is one of Lyapunov's 'critical' cases [133], [134] [132]. We apply 
a linear transformation to split the system into critical and non-critical variables; 
where the critical variables are those eigenvectors with zero eigenvalue and the non- 
critical variable are the others. Then we apply a non-linear transformation which 
will eliminate the influence of the critical variables upon the non-critical ones at the 
leading order. This gives {w, v) — > {W, W + V), and the system (2.4.24) becomes: 

W' = -2W{W + V) (2.4.25) 
V = -V + 2W{W + V) 

with the critical point at (M/, V) = (0,0). The Lyapunov procedure for the system 
(2.4.25) is to set the hnearly stable variable, V, equal to zero in the W equation so, 

W' = -2W^ 

and so the second-order analysis shows that critical point {W, V) = (0, 0) is unstable. 
In fact this critical point is a saddle point as can be seen from the numerical-phase- 
plane plot figure 2.5. The unstable part correspond to a non-physical range of the 
cosmological variables, since it gives w < 0. The stable part corresponds to the 
range of values where w > 0. With respect to our (approximate) exact solution, the 
unphysical case will correspond to the range of xq < 0, since they lead to a < 0. 
Hence the asymptotic solution (2.4.23) is the stable late-time behaviour of the dust 
universes with small (. The phase plane structure and the evolution of ip vs. \nt is 
shown in figure 2.5. 

2.5 Some General Asymptotic Features 
2.5.1 Models with asymptotically constant ip and a 

We have seen that ijj, and hence the fine structure constant, a, tends to a constant 
at late times in accelerating universes with power-law and exponential increase of 
the scale factor. We can establish a useful general criterion for this asymptotic 
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Figure 2.5: Numerical plots of the phase space in the {w,v) coordinates, and the ip 
evolution with x = logit) for n—^. The triangle is a saddle point for w < and a 
stable node for w > 0. 



behaviour to occur for general a{t). Suppose that as t ^ oo both sides of equation 
(2.2.2) tend to a constant (which may be equal to zero). Thus 

, f tdt f dt ^ ^ 

as t — > oo if 



a 

Then for consistency we also require, as t — > cxd, that 



B I ^ ^ (2.5.2) 



A^exp[-2C] ^A 



so A cannot be chosen to be zero. Thus in all cosmological models for which (2.5.1) 
and (2.5.2) hold there will be a self- consistent asymptotic solution as t — > oo of the 



2.6 General Analysis of the Phase Plane Bifurcations 



54 



form 



where B,C are positive constants. We see that the n > 2/3 and de Sitter cases 
satisfy the conditions (2.5.1) and (2.5.2) hence i/j is asymptotically constant. The 
dust (n = 2/3) case, for which ip —>■ oo as t ^ oo satisfies (2.5.2) but violates (2.5.1). 



2.6 General Analysis of the Phase Plane Bifur- 
cations 

In previous sections we have analysed the ip evolution equation (2.2.2) for a range 
of variables which are physically realistic and correspond to expanding universes. 
We will now analyse the whole range for variables of the system (2.4.17). As before 
we see there are two critical points in the {w,v) plane, at (0, —1 + ^) and ((1 — 
3n){^ — 1),0). Linearising (2.4.17) about {wci,VcJ = (0,-1 + ^) and (^02,^02) — 
((1 — 3n)(^ — 1), 0) we obtain the following characteristics matrices: 

/2-3n \ io (3n-l)(2-3„)\ 

y 1 1 - 3n y \l 1 - 3n J 

The characteristic matrices are non singular except when n — ^ or n — ^. In 
the non-singular cases the critical points will be simple and the system defined by 
these differential equations is structurally stable [130], and there will be no 'strange' 
chaotic behaviour outside the neighbourhood of the critical points. Hence, the 
linearised system will have the same phase portrait as non-linearised one in the 
neighbourhood of the critical points. 

The evolution, with respect to changing n, of the signs of the determinant and 
the trace of these two matrices is given in the table 2.1. This show us that there are 
always two critical points in our system, an unstable saddle and an attractor (which 
changes from a spiral to a node). 

The cases n = | or n = |, where the determinant of the characteristic matrices 
vanishes, lead to a bifurcation of codimension 1, in particular, of Saddle- Node type 
[131], since they correspond to points where the determinants of the characteristic 
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n 


Critical Points {wc, Vc) 






((3n- 1) (1 - ^) ,0) 


(0,1 - 4?) 
V ' 2 / 




f— oo; k) 


Saddle Point (non-physiccil) 


Unstable Node 




det Ml < 0, Tr Ml > 


det M2 > 0, Tr M2 


> 


1 
3 


Origin 


Axis 






det Ml = 0, Tr Ml = 


det M2 = 0, Tr M2 


> 


(-•-) 


Sta.blp Sniral 


Saddle Point 




det Ml > 0, Tr Ml < 


det M2 < 0, Tr M2 


> 


1 
2 


Stable Spiral 


Saddle Point 






det Ml > 0, Tr Ml < 


det M2 < 0, Tr M2 


= 


(L-i) 


Stable Spiral 


Saddle Point 




(let Ml > 0, Tr Mi < 


dot M2 < 0. Tr M2 


< 


3 
5 


Stable Spiral (node) 


Saddle Point 






det Ml > 0, Tr Ml < 


det M2 < 0, Tr M2 


< 


V5' 3/ 


Stable Node 


Saddle Point 






det Ml > 0, Tr Ml < 


det M2 < 0, Tr M2 


< 


2 

3 


Stable Axis 


Stable Axis 






det Ml = 0, Tr Ml < 


det M2 = 0, Tr M2 


< 


(|;oo) 


Saddle Point (non-physical) 


Stable Node 






det Ml < 0, Tr Ml < 


det M2 > 0, Tr M2 


< 



Table 2.1: Stability of all the critical points of the equation of motion for ijj. 

matrices change sign, det Mi = detM2 = 0. At these values of n the nature of the 
system will change. Cosmologically, these points represent a change in the behaviour 
of the time evolution of the fine structure constant as can be seen from the figures: 
2.6, 2.4, 2.1, 2.2, 2.5 , 2.3, which display the time evolution of ip. When n starts to 
grows from — cxo to | the two critical points slowly converge at n = |. For example, 
in the n = case where we may without loss of generality set a = 1, equation (2.2.2) 
has the exact solution 

a = ex.p[2ilj] = A-^ cosh^[AN^/'^{t + to)] 

where A, to are constants. This is an unrealistically rapid growth asymptotically, 
i/j Gct, caused by the absence of the inhibiting effect of the cosmological expansion. 
The case of n = 1/4 is shown in figure 2.6, which shows the phase space trajectories 
and the evolution of ip vs. Int. 
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Figure 2.6: Numerical plots of the phase space in the {w,v) coordinates, and the ip 
evolution with x — log{t) for n — \. The sign is a saddle point and the square 
is an unstable source. 

At n = I we are in the situation where the two critical points collapse into 
a unique one at the origin, creating a saddle-node bifurcation and a concomitant 
change in the behaviour and evolution of ip. As n keeps growing the single critical 
point splits into two critical points again. They move apart until the radiation value 
is reached, n = ]^ (Tr M2 = 0)- In this case, ipis a, asymptotically monotonic growing 
function of time, with some small oscillations near the Planck epoch. However, note 
that in our universe the asymptote giving an increase of if) behaviour with time is 
never reached before the dust-dominated evolution takes over [107], [108], [109]. 

As the universe evolves to the dust-dominated epoch, and n approaches the 
intermediate behaviour n — the two critical points start to coalesce again into a 
single point. When n = | is reached, becomes a strictly monotonically growing 
function of time. When n reaches the value corresponding to a dust-dominated 
universe, n — ^, another saddle-node bifurcation occurs. The two critical points 
collapse into a single one. Again there will be change in the behaviour of ip for 
larger values of n. Accordingly, when n > 2/3, the two critical points reappear once 
again and ip becomes asymptotically constant in value. 

Notice, that although a bifurcation is something that 'spoils' the smooth be- 
haviour of a system, in our case, that won't happen, due to the physical constraints 
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of our variables. In reality due to those constraints, the physical system will never 
'feel' the abrupt change at n = | and n — ^. This is also due to the fact that the 
attracting critical point always lies in the physical range of the variables, while the 
unstable one disappears form the physical system when the bifurcations occur, as 
can be seen from the phase plane plots. 



CHAPTER 3 



Gauge-Invariant Linear Perturbations of 
Varying- Alpha Cosmologies 



Sera este espelho a porta de entrada para a casa da verdade 
em que ainda ndo fomos convidados? 

- Ernesto Mota - 



3.1 Introduction 

In the previous chapter, using a phase space analysis of the equation of motion of 
the field t/j (1.7.8), we have shown that BSBM theories have a number of appeahng 
properties. They predict that, in a homogeneous and isotropic universe, there should 
be no variation of a during the radiation era and none during any present or late-time 
curvature or cosmological constant dominated era. During the dust era a should 
grow (leading to Aa/a < as observed) but only as log(i). Those results were also 
obtained by Barrow, Sandvik, and Magueijo [107, 108, 109, 110] using a different 
approach. 

This behaviour allows the quasar data to be accommodated without producing 
conflict with recent geonuclear limits on allowed variations of a, like the Oklo natural 



58 



3.1 Introduction 



59 



reactor limits of Aa/a < 10~^ [20, 12, 11, 22, 8] because they are imposed at a very 
low effective redshift of 2; 0.15, at which time the universe has begun accelerating 
and all variations of a are damped out. Recent deductions of a possible upper 
limit of Aa/a < 3 x 10~^ at 2; = 0.45 from nuclear (3 decays are potentially more 
restrictive [135, 25, 21]. These observations may be indeed problematic, since there 
is, at least, an order of magnitude difference between Aa/a from the quasar data 
and the /3-decay rate. 

A possible approach to solve the discrepancy of the two results is to claim that 
one of the measurements has a non-identified problem which was not taken into 
account in the observations. However, it must be remembered that both the nuclear 
limits are derived from a local solar-system environment. In the absence of a theory 
relating the value of, and rate of change of, a on the cosmological scales where quasar 
lines form to their values in the virialised local inhomogeneities where galaxies, stars 
and planets form, one should be wary of ignoring the possible corrections that must 
be introduced when comparing planetary and quasar bounds: the density of the 
Earth would not, for example, be a good indicator of the density of the universe. 
Adding to that, we know that any varying-a theory implies the existence of a field, 
responsible for variations of a, which is coupled to the matter fields. It is then 
natural, that strong variations of the density of matter in the universe, lead to 
spatial variations of ip. Thus, inhomogeneity is an important factor in the study of 
varying-constant cosmological theories. 

The first approach to address this problem, and the possibility of variations of a 
be infiuenced by matter inhomogeneities, is to consider the so called linear regime of 
the cosmological perturbations. In this regime, all the perturbations in the matter 
fields are very small. Much smaller than the background value of the fields. 

In this chapter, we are then going to study the evolution of small, gauge-invariant 
perturbations to the exact Friedmann-Robertson- Walker solutions of the BSBM 
theory. The results of this investigation will reveal whether one must worry about 
fast growth of small initial inhomogeneities in the value of a, which would lead to 
spatial variations of the value of the fine structure constant that might be more 
significant than the time variations at late times. 

Our discussion is organised as follows. In section 2 of the we give the gauge- 



3.2 Linear Theory of Cosmological Perturbations 



60 



invariant perturbation equations following the development used in general rela- 
tivity. In section 3 we specify the BSBM cosmology with varying a and derive 
the gauge invariant linear perturbation equations which couple the perturbations in 
the gravitational field to those in a and the density. In section 4 we solve for the 
time-evolution of small inhomogeneities in the fine structure constant on large and 
small scales for radiation, dust, and cosmological constant dominated expansion of 
the background universe. In section 5 we extend these studies into the non-linear 
regime by means of numerical solutions for flat and closed FRW universes. The 
evolution of spherical curvature inhomogeneities in density and in a is followed by 
computing the difference in time evolution between the FRW models of different 
curvature. These studies also reveal for the first time the behaviour of a in closed 
FRW models in the BSBM theory. 

Units will be used in which c = h = 1; Greek indices run form to 3 and Latin 
indices only over the spatial degrees of freedom 1 to 3. The Einstein summation con- 
vention is assumed; a{t) is the scale factor of the background Friedmann- Robertson- 
Walker (FRW) metric and G is Newton's gravitation constant. Recall that in the 
BSBM varying-a theory, the quantities c and h are taken to be constant, while e 
varies as a function of a real scalar field ip, with a — e^, hence 
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Observations indicate that the universe is nearly homogeneous and isotropic on 
large scales, in particular at early times where metric perturbations to the cosmic 
microwave background radiation are observed to be small ~ 10~^ << 1. It is usually 
assumed that there exist small primordial perturbations which slowly increase in 
amplitude due to gravitational instability to form the structures we observe at the 
present time on the scales of galaxies and clusters of galaxies. The problem of 
describing the growth of small perturbations can be tackled using the cosmological 
perturbation theory, which is based on expanding the Einstein equations to linear 



a 



e 



eoexp[V'] 
cco exp[2'0] 



(3.1.1) 
(3.1.2) 
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order about the background metric. The theory was initially developed by Lifshitz 
[136], but many other approaches and improvements were done until today, see for 
instance [137, 138, 139, 140] for reviews. 

Here we will follow the gauge invariant formalism developed in [141, 142] and we 
refer the reader to this references for more details in this particular approach. Also, 
we will a priori assume the existence of small inhomogeneities at some initial time 
in the early universe. Hence we will not discuss any particular model to create the 
primordial fluctuations in the matter fields, see [143, 144, 145, 146, 147] for different 
proposed models and hypothesis. 

3.2.1 The background and metric perturbations 

Assuming that the metric of spacetime deviates only by a small amount from a 
homogeneous, isotropic FRW spacetime, which is defined to be the background 
universe. It is convenient to split the metric into two parts: the first is the back- 
ground metric, which represents the ideal background spacetime, and the second 
is its perturbation, which describes how the 'real' spacetime deviates from the ide- 
alised background. This approach is mathematically well defined since it was shown 
in [148] that in the case of a FRW universe, solutions of the linearised field equations 
can be viewed as linearisations of solutions of the full non-linear ones. 
The background line element is: 



ds^ =W Qi^^dx^'dx" = - o?{t)'-iijdx^dx^ = a^{r]){drf' - 'yijdx'dx^), (3.2.1) 



where K = 0,1,-1 depending on whether the three-dimensional hypersurfaces of 
constant r) are flat, closed or open. The Einstein equations are: 



where rj is the conformal time. 



drj — a ^dt. 



We choose the background metric to be the FRW metric, so 




(3.2.2) 
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where is the Ricci tensor, R = R'^is the Ricci scalar and is the total energy- 
momentum tensor. We shall exploit the fact that the BSBM theory for varying 
a can be expressed as general relativity with a particular linear combination of 
energy-momentum tensors. For the moment consider the presence of a single energy- 
momentum tensor. In the next section this will be decomposed appropriately. 

For the background metric in equation (3.2.1), in conformal time, the Einstein 
equations reduce to the — equation 

7^2 ^ ^^2^0 _ (3 2.3) 

where 7i = ^ using conformal time, and the i — i equation: 

n' + = ^^a^T -K, T = TI^ (3.2.4) 

For the background metric (3.2.1), the space-space part of the Ricci tensor 
i?*- is proportional to 5j. Thus, for an isotropic background universe, the energy- 
momentum tensor must also be spatially diagonal, Tj oc 5], in order that the Einstein 
equations are satisfied. Differentiating (3.2.3) with respect to 77 and subtracting 2a' 
we get the continuity equation for matter VjuT(f = : 

dT^ = -(4T° - T)d\na. (3.2.5) 



We now introduce small perturbations around the FRW background and follow 
the gauge-invariant approach of Mukhanov et al [142]. The full line element may 
be expressed as: 

ds^ = ^^'^g^^dxi'dx" + Sg^^dx^dx", (3.2.6) 

where Sg^i, describes the perturbation. The full metric has been decoupled into its 
background parts and perturbation parts: 

g^u = + (3.2.7) 

The metric tensor has 10 independent components in 4 dimensions. For lin- 
ear perturbations, the metric perturbation can be divided into three distinct types: 
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scalar, vector and tensor, according to their transformation properties on spatial 
hypersurfaces [141, 149]. The reason for splitting the metric perturbation into these 
three types is that they are decoupled in the linear perturbations equations and 
so evolve independently. Neither of the vector and tensor perturbations exhibit 
growing instabilities in dust and radiation universes . Vector perturbations de- 
cay kinematically in an expanding universe whereas tensor perturbations lead to 
gravitational waves that do not couple to the isotropic energy-density and pressure 
inhomogeneities [142]. However, scalar perturbations may lead to growing inhomo- 
geneities which, in turn, have an important effect on the dynamics of matter and 
thereby on the time and space variations of the fine structure constant in the BSBM 
theory. In this thesis we will then consider only the scalar perturbation modes. 

It is shown in the literature, see for instance [141] , that scalar perturbations can 
always be constructed from a scalar quantity, or its derivatives, and any background 
quantity, such as the 3- metric 'jij. We can then construct any first-order metric 
perturbation in terms of four scalars 0, x, B and E, which are functions of space 
and time coordinates [142]: 



where \i represents the three-dimensional covariant derivative. The most general 
form of the scalar metric perturbations is then constructed using the four scalar 
quantities^, x, B and E: 



From the above equation and equation (3.2.7), the line element for the back- 
ground and for the scalar metric perturbations is 



ds'^ = a\r)){{l + 2^)drf - 2B\idx'dr) - [(1 - 2x)-fij + 2E\ij]dx'dx^)}. (3.2.9) 



-2a^(t), 
5gio = a^B^i, 
-2a^ {xiij - E\ij) 
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3.2.2 Gauge-invariant variables 

The homogeneity of a FRW spacetime gives a natural choice of coordinates. But 
in the presence of first order perturbations we are free to make a first order change 
in the coordinates, that is, a gauge-transformation. This gauge freedom leads to 
several complications. For example, not all perturbed metrics correspond to per- 
turbed spacetimes. It is possible to obtain an inhomogeneous form for the metric 
in a homogeneous and isotropic spacetime, by a choice of coordinates. Hence, it 
is important to be able to distinguish between physical inhomogeneities and mere 
coordinate artifacts. 

The gauge problem of the cosmological perturbations is intimately related to 
the fact that we actually deal with two spacetimes, the physical space time, corre- 
sponding to the real universe and a fictitious one, which defines the unperturbed 
background. We need then to establish an one-to-one correspondence between the 
two spacetimes. For instance, when a coordinate system is introduced in the back- 
ground we need to have a correspondent coordinate system in the real universe and 
vice- versa. This is necessary in order for us to be able to define a background space- 
time into the real universe and to compare quantities between them. The need for 
that is clear from the definition of a perturbation in a given spacetime: a pertur- 
bation in a given quantity is the difference between its value at some event in the 
real spacetime and its value at the corresponding event in the background. This 
means that even if a quantity behaves as a scalar under coordinate transformations, 
its perturbation will not be invariant under gauge transformations. As a result we 
might end up with spurious gauge modes in the solutions, which have no meaning 
whatsoever. 

There are two essential ways to deal with the gauge problem. The first is to 
choose a particular gauge and compute all the quantities in that gauge. If the 
gauge choice is well motivated, the perturbed variables will be easy to interpret. 
However, the task of selecting the best gauge for any given situation is not always 
trivial. See for example [136, 52] for more details on this approach. The second is 
to describe perturbations using gauge-invariant variables. We will follow here the 
second approach. 

Gauge-invariant variables are unchanged under all infinitesimal scalar coordi- 
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nate transformations, so they are independent of the background coordinates. Such 
quantities can be constructed out of the four scalar functions E and B, see 

[141]. The simplest gauge- invariant linear combinations which span the space of 
gauge-invariant variables that can be constructed from the metric variables alone 
are [142]: 

^^^^[{B-pal^ ^ = « (3.2.10) 

Note that any combination of gauge-invariant variables will also be gauge invari- 
ant. In general, a scalar quantity /(?], x) defined in the spacetime can be split into 
its background value and a perturbation f{r], ) + Sf{ri, x). Since, in general, Sf is not 
gauge invariant, we cannot use this scalar quantity without modification if we want 
to have gauge-invariant equations. Hence, we consider the following gauge-invariant 
combination: 

5f<^')^5f + f,{B-E'). (3.2.11) 



The freedom of gauge choice can be used to impose two conditions on the four 
scalar functions. The longitudinal gauge is defined by the conditions E — B — 0. 
This gauge choice has the advantage of ruling out the complications of residual gauge 
modes. Also, in this gauge and x coincide with the gauge-invariant variables $ 
and respectively. In this longitudinal gauge, the metric takes the form: 

ds^ = a\r))[{l + 2$)(i?7^ - (1 - 2^If)^ijdx'dx^)], 

and the gauge invariants $ and ^ become the amplitudes of the metric perturbations 
in the longitudinal coordinate system. In the case where there are no space-space 
components in the energy-momentum tensor, so Tj oc we have that $ = \I' and 
there remains only one free metric perturbation variable which is a generalisation 
of the Newtonian gravitational potential. As can be seen from the above equation, 
the gauge invariants $ and ^ have a simple physical interpretation: they are the 
amplitudes of the metric perturbations in the longitudinal coordinate system. 
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3.2.3 The linearly perturbed Einstein equations 

For small perturbations of the metric, the Einstein tensor can be written as = 
^°-'G(^+(5G'(^, and the energy-momentum tensor can be split in a similar way. The fully 
perturbed Einstein equations can be obtained for scalar perturbation modes with 
the line element given by (3.2.9). However, these equations are not gauge invariant, 
since they contain non gauge-invariant quantities. In order to have gauge-invariant 
equations we need to replace (f) and x by the gauge-invariant variables $, ^' and 
B — E', and to construct the gauge-invariant equivalents of and 5Tj^, we need 
to rewrite then as [142]: 

SG^f = 5(7} + ((°)G'}y(S -£;'), (3.2.12) 
5Gf^' = ^ ((0)^0 - 1 (o)G^)(5-^')|., 

o 

and analogously for ^T^, 



^^^gi^o ^ ^^o^^(o)^oy(^_^/)^ (3 2.13) 

^ 5T; + {^''^T;y{B-E'), (3.2.14) 



^^(gi)o ^ 57;o + ((o)tO-1 (o)rfe^)(5-E')|i. (3.2.15) 

3 



The components of the perturbed Einstein equations linearised around small 
perturbations of the background are SGI^^^ ^ — SnGST^i^^^ 

= - 3n^' - 3^n^ -k)= ATrGa^ST^''^'' 

= a,(a$)' = 47rGa5T/^')° (3.2.16) 
SGi = ^" + 3W + ^{2n' + n^-K)^-4TTGa^STj''^' 

where we have already simplified the equations since the energy-momentum tensor 
that we will be considering in section 3 has no space-space components and so 

In order to close our system of equations, we need equations of motion for the 
matter formulated in a gauge-invariant way. This requires explicit expressions for 
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the energy-momentum tensor and so we must now specify the BSBM theory. 

3.3 Cosmological Perturbations and the BSBM 
Model 

It was shown in chapter 1 that, in the context of the BSBM model, the homogeneous 
evolution of ip does not create significant metric perturbations at late times and the 
cosmological time-evolution of the expansion scale-factor is very well approximated 
by the usual power-laws found in k, — FRW models filled with a perfect fluid (see 
also[110] ). 

Therefore, we will assume there are no major modifications in the perturbed 
spacetime which would lead to changes of the behaviour of the perturbed variables 
of a perturbed FRW spacetime filled with perfect fluid. That is, we will assume that 
the energy-density perturbations and the metric potential are the same as in a FRW 
universe with no variation of a. Physically, this is to be expected for most of the 
evolution, although this assumption might break down (along with much else) on 
approach to initial and final cosmological singularities. It is a reflection of the fact 
that the changes in a have negligible feedback into the changes in the expansion, 
which arc governed to leading order by gravity. The principal effects arc those of 
perturbations in the matter density and expansion rate on the evolution of a. This 
simplification will allow us to write the time and space variations of the scalar field. 
Sip, as a functions of pm, Pr, Spm, Spr , $ and ip, where Spm, Spr and $ will be given 
by the solutions found in ref. [142] for universes with no variation of a; the field 
ip will be given by the solutions found previously in chapter 2. In order to find an 
expression for Sip in terms of these quantities we need to write the gauge-invariant 
linearly perturbed Einstein equations for the BSBM model. 

3.3.1 The background equations 

We vary the action (1.7.2) with respect to the metric to obtain the generalised 
Einstein equations: 
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where — y= \^i_,J^ - is the energy-momentum tensor for perfect- fluid matter 
fields, and 

where -u^ = 5q is the comoving fluid 4-velocity; ^ and T^™ ^ are the energy- 
momentum tensors for the kinetic energy of the field ip and the electromagnetic 
field respectively: 

Note, the total energy density of the electromagnetic field is the sum of the Coulomb 
energy density C^p„i and the radiation energy density pr , where — 1 < C ^ 1 i^ the 
fraction of mass density of matter in the form of the Coulomb energy. We will 
then consider T^'^ as a perfect fluid: 

= i\C\Pm + Pr+Pm+ Pr)e~'^''' U^U^ - {pm + Pr) C^^'^fi'^v 

The propagation equation for ^0 comes from the variational principle as: 

d, [V^g^'^d^^l^] = -i^e-2^£em. (3.3.1) 

to 

This equation determines how e, and hence a, varies with time. 
The background equations can now be explicitly obtained: 

?,7f = SttGo^ (^p^ + (p, + |C|pJe-2^ + |a-V') +«'A-3«; (3.3.2) 
m' = -AnGa" (pm (1 + |C|e"''^) + 2p,e-'^ + 2u;a-V') + a'K 
where A is the cosmological constant. The equation of motion for the '4> field is: 

27111^' + il^" = ?^p^e-^^ (3.3.3) 

U) 

The conservation equations for the matter fields, pr and pm respectively, are: 

p'r^ + mpm = (3.3.4) 

p'^ + mp'^ = 2^'pr. (3.3.5) 
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3.3.2 Linear perturbations in the BSBM theory 

The perturbed components of the total energy-momentum (T^°*"' = T^"* + T^^, + 
T^™) arise from perturbations of the different matter fields which are time and space 
dependent. In particular, we have Pm Pm + Spm, Pr Pr + ^Pr and ip ^ ip + Sip. 
Note also that we have to perturb the fluid 4- velocity field, so we have Ui — > Ui + 6ui, 
where i = 1,2, 3. 

In order to have gauge-invariant equations we need to express the perturbed 
energy-momentum tensor in terms of the gauge-invariant energy density, pressure, 
scalar field and velocity field perturbations . The gauge invariants 5pm^\ 5p^r^^ and 
are defined in the same way as the gauge-invariant perturbation of a general 
four-scalar, see equation (3.2.10), so: 

5p^^^ = 5pm + p'^{B-E'), 5p^^'^ = 5p + p\B-E'), 5i^^^'^ = 5^; + i;\B - E'), 

where 5p is the perturbed pressure for a specific component and the gauge-invariant 
three- velocity 5u^^^^ is given by [142]: 

5uf^ ^5ui + a{B-E%. 

The physical meaning of the quantities which enter the gauge-invariant equations is 
very simple: they coincide with the corresponding perturbations in the longitudinal 
gauge. From now on we will drop the superscript {gi) since we will always be dealing 
only with gauge-invariant quantities. 

We can now write the gauge-invariant linearly perturbed energy-momentum ten- 
sor: 

5pm, (3.3.6) 
{pm+Pm)a~^5ui, (3.3.7) 
-5pn,S], (3.3.8) 



^rpmatter 



5T- 



matter 



5T, 



matter i 
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STf" = a;a-2 (V'^V^' - , (3.3.9) 
STf" = ua-^i)'dM, (3.3.10) 
5Tf' = a;a-2 (^/$ - VV) (^j, (3.3.11) 



3 



5Tr' = i\C\5pm + 5pr)e-'^-2e-^^{\C\pm + Pr)5iP, (3.3.12) 
STr"" = e-^'^{\C\prr, + pr+p„,+Pr)a-'6u.„ (3.3.13) 
^j:emi ^ ^26ilj{pm+Pr)e-'^-{6pm + 6pr)e-^^)S}. (3.3.14) 

We have assumed there are no anisotropic stresses in the energy-momentum ten- 
sors and we have considered only adiabatic perturbations; that is, we consider the 
pressure perturbations to depend only on the energy-density perturbations. In the 
dust and radiation cases this means p^ — 0, Sp^ — 0, or Pr — \pr and 5pr — \pr-i 
respectively. 

The fully perturbed gauge-invariant Einstein equations can be obtained from 
(3.2.16) using the expressions above for ( 3.3.6-3.3.14). We have 

$ (^-37^^ + 3A; + 4G'7ru;V'^) + - AGt^ujiI)' d^)' - SW = 

iGna'^e-^^ [(e^^ + |C|) + Spr - 25^ (p, + |C|p™)] (3.3.15) 

7^V'$ + V'$' = 4G'7ra;V'v'(^V' - ^ae-'^ [4p, + 3 {e^^ + |CI) Pm] V5i^ 

$ (n^ + 2n' -k + 4G7ra;^'^) + 3W + $" = 

^ [a'e-'^ {Spr - 2prS7P) + SujtP'StP'] (3.3.16) 

It is useful to write the perturbed energy-momentum conservation equations for 
each component: 

Sip" = ^e^^^a^ [6pm + 2pm ($ - 5ip)] - 2nSip' + 7^Sip + iip'<^' , (3.3.17) 
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Sp'^ = ~ 6pr {en - 3V^') + Pr (2v6u - 36^)' - 6$') , (3.3.18) 
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5p'^ = -3n6p^ - p„ (v5u - 3$') . (3.3.19) 

From these expressions it is clear that perturbations in a are sourccd by pertur- 
bations in the dust, but not vice versa. Hence we expect that, in a dust-dominated 
universe with varying a, the cold dark matter perturbations will behave as in a 
dust-dominated universe with no varying a. Notice however that the same cannot 
be concluded so easily for a radiation-dominated era, since there is a source term in 
equation (3.3.18) proportional to Sip' . We expect this term to be negligible at large 
scales, but that might not be the case on small scales. 

The gauge-invariant perturbation for a is given by equation (3.1.1) as 

a 

From equations (3.3.15), (3.3.16), (3.3.16), using (3.3.2) to simplify, we obtain the 
general form for 6ip, the perturbation to the scalar field which drives variations of 
the fine structure constant, as 

{5pm - 2p^$)] + 'ie^^ 2$ (m'^ -k- AGtiujiP'^'^ - V^* + m^' + ) 

This is a gauge-invariant expression for dip, written as a function of the gauge- 
invariant quantities $, 5pm-i and Spr- 

3.4 Evolution of the Perturbations 

In a spatially flat universe, k = 0, the general gauge-invariant expression for Sip 
becomes: 

= 8Grf(2p!+3|C|p.,.) (^''^"^ - ^ ' ^ 1^1' 

{5pm - 2p^$)] + 362^^ $ - SGtto;^'^) - V^^ + m^' + ) . 

It was shown in chapter 2 that the cosmological evolution of the metric scale 
factor is unchanged (up to very small logarithmic corrections) to leading order by 
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the time evolution of ijj. The dominant effect is that of the evolution of the scale 
factor on the evolution of ip through its propagation equation. 

It is known that perturbations of massless scalar fields, or scalar fields with a 
very small mass are negligible with respect to perturbations in the matter fields 
and the gravitational potential. Guided by this, in order to obtain the evolutionary 
behaviour for we will assume that the matter field perturbations, Spm and Spr, 
and the metric perturbations, are unaffected by the perturbations to leading 
order. We will assume that these three quantities will therefore be the same to 
this order as they are in a fiat FRW universe filled with barotropic matter and a 
minimally coupled scalar field. These assumptions are valid if the energy density of 
■0 is much smaller than the energy density of the matter fields, so $ will be driven 
only by the matter perturbations. If we examine the perturbations in the non-linear 
regime it is confirmed that ^ ^ ^ (see figure (3.4)). 

3.4.1 Radiation-dominated universes 

In any radiation-dominated era in which the expansion of the universe is dominated 
by relativistic particles with an equation of state p — ^p, we can neglect the non- 
relativistic stresses in the universe, in particular, the cold dark matter and the 
cosmological constant since pr » Pm » Pa to a good approximation. If we 
assume that Pa — Pa — — k and 5pm — 0, the background equations of motion 
give the usual conformal time evolution for the scale factor and the energy density 
of the radiation: 

Pr — ProCT^^^^ Cb — s/ '&'^G ProT] 

where p^o is a constant. 

The perturbations in the barotropic matter fiuid and the potential $ come from 
the equations (3.3.15), (3.3.16). The usual fiat FRW solutions with constant a 
are obtained by setting the terms proportional to to zero, (3.3.16). The general 
solution of these equations can be obtained by expanding the physical quantities in 
terms of the eigenfunctions of the operator (where —k"^ denotes the eigenvalue 
of this operator) and solving for each mode separately. Resuming the terms, the 
general solutions of the linearly perturbed equations for the potential $ and the 
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barotropic matter are: 

$ = r]~^{[wr]COs{wr]) — sm{wr])] Ci + [wrj sm{wrj) + cos{wr])] C2}e^^ (3.4.2) 

and 



Pr V 



1 - ^ {7]wf 



cos{r]w) + [(jyw)^ — l] 8in{r]w)} (3.4.3) 



+{ [l — (rjw)'^] cos{r]w) + rjw 



1 - \ ivwf 



sm{r]w)}C2)e 



jkx 



where we have expanded the general solution in plane waves since we are assuming a 
spatially flat universe; Ci and C2 are arbitrary functions of the spatial coordinates; k 
is the wave vector mode and w = k/ a/3- Note that these quantities are all expressed 
in their gauge-invariant format. Finally, to calculate the explicit time dependence 
of Sip, we need to use the background solution for -0: 

V' = ^log(8iV) + llog(|r^2). (3.4.4) 

This was found in chapter 2 and [109, 1] for a radiation-dominated universe, where 
N — —^pjnCL^ > since C < in the magnetic energy dominated theories considered 
by BSBM. It is important to notice that in universes with an entropy to baryon ratio 
{S ~ 10^) hke our own, t/j does not experience any growth in time [108]. The constant 
term on the right-hand side of (3.4.4) dominates the solution for ■0(77) throughout 
the radiation era. Numerical solutions confirm this freezing in of ■0, and hence of a, 
during the radiation era, [107]. 

Large-scale perturbations in a radiation-dominated era 

In the long- wavelength limit (^77 << 1) where the scale of the perturbation exceeds 
the Hubble radius, we can neglect spatial gradients. So, in this limit, ^ becomes: 

S^^-—cx — e^'^^V^Cs - 47rG'we^^^C27?-^ 
2 a 2?7 

From this expression we can sec at that on large scales the inhomogcneous per- 
turbations in a will decrease as a power-law in time. This behaviour agrees with 
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the one found in [110] by other methods. 

Small-scale perturbations in a radiation-dominated era 

On scales smaller than the Hubble radius {wrj » 1) the dominant terms are pro- 
portional to V^Ci and V^C2, so the asymptotic behaviour for ^ will be: 

1 Sex 1 

Sijj = - — oc --W [cosfryw) V^Ci + sinfryw) V^Ca] e^"" + (3.4.5) 

2 a 2 

+^ rsin(7^w)V^Ci - cos(ryw)V^C2] e^^"" 
2ri 

On small scales we can see the perturbations on a will be oscillatory. This 
behaviour is new and does not coincide with the ones found in [110]. 

In reality, wc expect dissipation of the adiabatic fluctuations to occur by Silk 
damping [150] and small-scale fluctuations in a will also undergo decay as their 
driving terms damp out. However, while the exact solution for the evolution of a 
is a linear sum of a constant and a slow power-law growth, the power-law evolution 
does not become dominant by the end of the radiation era in universes like our own 
with entropy per baryon O(IO^). 

3.4.2 Dust-dominated universes 

In the case of a flat dust-dominated universe, flllcd with a p = fluid, we can assume 
Pr = Pa = = K and 5pr = 0. Again, in order to obtain an explicit expression of 
Sip from equation (3.4.1), we will assume that the mattcr-flcld perturbations, pm 
and 6pm, and the metric perturbation, $, are not affected by the ip perturbations to 
leading order. Thus, we will assume that these functions behave as in a perturbed 
flat FRW dust universe. 

In general, for a flat dust universe we have the following background solutions: 

Pm = PmoO , a = 77 , 

where is constant. 

As before, we can calculate the most general gauge-invariant solutions for the 
energy-momentum perturbations Spm, and for the potential $, assuming that the 
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field does not affect them to leading order, so their time dependences are [142] : 

$ = Ci + C2V~^ (3.4.6) 

and 

^ = ^ [(^'v^Ci - 12C,) + {v\'C, + 18C,) V-'] , (3.4.7) 

where Ci and C2 are arbitrary functions of the spatial coordinates. 

Once again, we need the background solution for ■0, in order to calculate Sip as 
an explicit function of the conformal time. We use the asymptotic solution 

V^=llog[27Vlog(^r;3)], 

which was found in [110], as an asymptotic approximation for a dust-dominated 
universe which is in good agreement with numerical solutions, and where, as above, 
N — — '^pjnO? > is a constant. 

Large-scale perturbations in a dust-dominated era 

On scales larger than the Hubble radius we can neglect the terms proportional 
to 9jCi,V^Ci, V^C2 and diC2] so in this hmit we have the following asymptotic 
behaviour for the non-decaying mode: 

iV, = 1^^^ oc -2C, - ^-^C, (3.4.8) 

Therefore, on large scales the inhomogeneous perturbations of a will not grow in 
time by gravitational instability. This behaviour can be understood with reference 
to the general evolution equation for ip in Fricdmann universes. On large scales, 
where spatial gradients in ip can be neglected with respect to its time derivatives, 
we may view inhomogeneities in density and in ip as if they are separate Fricdmann 
universes of non-zero curvature (/t 7^ 0). The growth of inhomogencity can be 
deduced by comparing the evolution of a in the /t 7^ universes with those in the 
K = model (a more detailed numerical study of this model will be presented in 
section 5 below). In effect, this uses the BirkhofF-Newton property of gravitational 
fields with spherical symmetry [151]. We note that the (3.3.3) evolution equation 
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has the simple property that ip cannot have a maximum because > when > 
because N > , [110]. This result holds irrespective of the value of k = 0. Thus 
ip and a will continue their slow increase in both over-densities, under-densities 
and the fiat background until we reach scales small enough for spatial derivative to 
come significantly into play. This has the important consequence that we do not 
expect large spatial inhomogeneities in a to have developed. However, it should be 
noted that the sensitivity of the observations of varying-a effects in quasar spectra 
is sufficient to discern variations of redshift space smaller than 0(10""^), which is of 
the same order as the amplitude of density fluctuation on very large scales in the 
universe. 

Small-scale perturbations in a dust-dominated era 

On scales smaller than the Hubble radius the dominant terms are those proportional 
to V^Ci and V^C2, so the asymptotic behaviour for the growing mode will be: 

6iP oc ^V^^^Ci (3.4.9) 

This shows that perturbations of a will grow on small scales. This result is a product 
of the assumption that on small scales the universe can be considered as being filled 
by an homogeneous and isotropic fluid, however we know this is not true below the 
scale where gravitational clustering becomes non-linear. On these small scales we 
also have to worry about new consequences of inhomogeneity which have not been 
included in our analysis. For example, the constant parameter N — —^^o,^ ^ C^m 
will vary in space due to inhomogeneity in the background matter density parameter 
Qjn and in the dark matter parameter —1 < C < 1- We have assumed that C < for 
the cold dark matter on large scales in order for the cosmological consequences of 
time-varying a to be a small perturbation to the standard cosmological dynamics. 
But on small scales the dark matter will be baryonic in nature and so C > there. 
Hence, we expect C to be signiflcantly scale dependent as we go to small scales. 
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3.4.3 Accelerated expansion 

In an era of accelerated expansion, a > 0, as would arise during inflation or during 
a A- or quintessence-dominated epoch at late times, we can consider the scale factor 
to evolve as a power-law of the conformal time as a = r^"", where n > 1 and r] runs 
from — oo to 0. The case of a = rj^^ corresponds to a A -dominated epoch. 

As in the previous sections, we will assume that all the other matter components 
which fill the universe will behave exactly as in a perturbed FRW universe with no 
variations of a. We also assume that neither of the dust and radiation perturbations 
will affect the behaviour and evolution of dip, and that these perturbations are 
negligible with respect to the A stress driving the expansion, so we will consider 
Spm = and 5pr — 0. 

We saw in chapter 2, that a universe which is undergoing accelerated expansion, 
the asymptotic solution for is a constant, so we will assume that ip — in the 
background, where ■0oo is a constant. Thus, equations (3.3.15) and (3.3.16) become: 

8Gn\C\6^Pm ^ ^ 3n (r^$ - n$') ^ ^ 



n 



(2 + n) $ + 77 (ji^" - 3n$') = 



where we have also considered pr — 0, but Pm ^ ^ because of the coupling with 
in the equation of motion of the scalar field. Note that if we had also set Pm — ^ 
here, we would have imposed a no a-variation condition: 5%l) — 0. 
Integrating the last equation, we obtain 

where Ci and C2 are arbitrary functions of the spatial coordinates. Note that since 
n > 1 in accelerating universes we see that $ will decay in time as 77 — >■ 0. From 
this solution for $, we obtain: 

15a 6^^^°° -3+n-Vl+n(-2+5n) V / . -\ 

= o— = -T7^ ^ 3n 1 + n- Vl + n 5n-2) Ci+ 



3n 



(1 + n + Vl + n(5n-2)) 77^1+^(5^-2) C2 + 2rfsj''Cr + 2ri'^^/^+^^^^V^C2 
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which is a decaying function of the conformal time when n > 1 and a constant when 
n — 1. Thus, in accord with the expectations of the cosmic no hair theorem, the 
universe approaches the FRW model and a is asymptotically constant at late times 
[152]. 

Large-scale perturbations during accelerated expansion 

On scales larger than the Hubble radius we can neglect the terms proportional to 
the spatial derivatives; so in this limit we have the following asymptotic behaviour: 

S^l^ oc -— —T] ^ 1 + n- Vl + n(5n-2) Ci 

16Gp„,„7r|C| LV / 

+ (l + n+^l + n{5n- 2)) 77 Vi+"(5«-2)C2 (3.4.10) 

Therefore, as expected, on large scales during an accelerated era inhomogeneities in 
a will decrease on time when n > 1 and will be a constant when n—1. 

Small-scale perturbations during accelerated expansion 

On scales smaller than the Hubble radius the dominant terms are the ones propor- 
tional to V^Ci and V^C2, so the asymptotic behaviour will be 

6^'*'^°° l+n-Vl+n(5n-2) V „ /T" 7; „ 1 , ^ 

5V'oc-— -'n V^Ci + TyVi+^Mv^Cs (3.4.11) 

This confirms that perturbations of a, as 77 — > 0, will decrease on small scales when 
n > 1 and will be constant when n—1. 



3.4.4 Summary of behaviour 



In Table 3.1 we summarise the time evolution of small inhomogeneities in a found 
under different conditions in this section. 

We should point out that those results were obtained under the assumption 
that there is no back-reaction, in the matter fields perturbations, due to the per- 
turbations in the fine structure constant. Back-reactions may have an important 
effect, and should be taken into considerations, specially, when there are coupled 
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Universal equation of state Time Evolution of the perturbations 



__L 2 ft 





Large scales 


Small scales 


Radiation- dominated epoch 


Decaying 


Oscillatory 


p = ^p, (I .3C // 






Dust-dominated epoch 


Constant 


Growing 


p = 0, a (xrf' 






Accelerated expn a oc ?7~", n > 1 






A-dominated 


Constant 


Constant 


p — —p, a oc r)''^ 






Power-law acceleration, n > 1 


Decaying 


Decaying 


p — wp, w < 0, a oc 77~" 







Table 3.1: Time evolution of small inhomogeneities in a found for the different 
epochs the universe went through. 



fluid-perturbations in the perturbed Einstein equations [153, 154, 155]. Perturba- 
tions in a given fluid, and so in its energy-momentum tensor, affect the dynamics 
of spacetime. These perturbations would then affect the evolution of the different 
components in the universe. Although, if the perturbations in a given fluid are very 
small and do not change much the spacetime geometry, then we can neglect its 
effects up to a reasonable approximation. 

Varying-a models have a scalar field coupled to the matter, a proper study of the 
possible back-reaction should then have been done. Although, there are very good 
indications that perturbations in a will not affect at all the matter inhomogeneities. 
The reason is the fact that the background universe is not affected by the inclusion 
of the scalar field responsible for the variations of a, because the energy density 
associated to the scalar field, when compared to the matter fields, is very small. 

Nevertheless, this assumption should be checked. In order to do that, we will 
use a simple toy model to mimic the situation when the perturbations in ip are large 
and so might affect the evolution of the other perturbations. We will study what 
is called the non-linear regime of the cosmological perturbations. In this regime, 
the approximation of using a first order perturbations, about the background, is 
no longer valid. Both the matter and perturbations will be no longer considered 
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small with respect to its background values. 

3.5 The Non-Linear Regime - A Toy Model 

In order to study the evolution of inhomogeneities in a beyond the domain of linear 
perturbation theory we need to use a different model. The simplest approach is to 
confine attention to spherically symmetric inhomogeneities. This will be done by 
comparing the solution of the BSBM theory for a in a closed {k — 1) universe, with 
the solution for a in a flat {k — 0) universe. We are assuming a Birkhoff property 
for the BSBM theory so that we can treat the perturbation as an independent closed 
universe [151]. This is a standard technique in general relativity which was flrst used 
by Lemaitre [156]. 

We define the alpha 'over-density perturbation' (which is not necessarily small) 

by 

6a ^ a^=i - a^=o 
ot q;k=o 

where is the solution of equation (3.3.1) for a universe with curvature k. 
3.5.1 Radiation-dominated era 

There is no non-linear structure formation during the radiation epoch. Nevertheless, 
it is interesting to see what would be the magnitude of the perturbations in a 
with respect to the matter fields. The scale factor for a radiation-dominated closed 
{k = 1) FRW universe is given by a = sin(77); for a fiat {k = 0) FRW universe the 
normalised scale factor is given hj a = rj. 

The evolution of a can be seen in figure 3.1 along with the evolution of a for 
fi: = 1,0. 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Figure 3.1: The evolution of a{r]) and a{r]) for radiation- dominated universes with 
K — (dashed) and k — 1 (solid). 

As expected, we can see there is no difference in the evolution of a at early times 
and q; oc 77 as it was found in [1]. When the difference between the scale factors of 
the two universes becomes significant, the behaviour of q;«;=i begins to deviate from 
that of aK=o- We see that clearly grows faster than aK=o- The difference in 
the growth rates become very significant near the expansion maximum of the bound 
region (rj n). However, after this time our assumption that the background is not 
affected by changes of ip in the cosmological equations that describe the background 
universe breaks down, since the kinetic energy of the scalar field will diverge and 
can no longer be neglected in the Friedmann equation. We expect the behaviour 
near the final singularity to be similar to the kinetic-dominated evolution near the 
initial singularity discussed in ref. [108]. 

From figure 3.2 we see that the variations of a will become increasingly important 
as f] approaches the second singularity. Notice that although there is a considerable 
growth in the perturbations in a, when we compare them with the perturbations in 
the radiation, they are not as significant as can be observed from the evolution of 
the ratio ^ versus 77 in figure 3.2. 
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Figure 3.2: The evolution of ^/-^ vs rj in radiation- dominated universes. 
3.5.2 Dust-dominated era 

During the dust-dominated phase of a closed universe [k = 1) the normahsed scale 
factor is given by a = 2 (1 — 005(77)), while for a fiat universe {k — 0) the normalised 
scale factor is given hy a — rf. Integrating equation (3.3.1) for both cases we obtain 
the evolution of a for both cases. 

The evolution of a{r]) can be seen in figure 3.3 along with the evolution of a{r]) 
for K = 1,0. As in the radiation case, we can see there is little difference in the 
evolution of a at early times, since the scale factor for the closed model evolves very 
similarly to the flat one for 77 << 1 and a oc ln(77). The differences between a^=i 
and q;«;=o cases start to appear when 77 ~ 1, when the nonlinear regime commences. 
These differences become more accentuated near the second singularity, but once 
again this is the region where our approximations break down. 

Notice that although there is a considerable growth in the perturbations in a, 
when we compare them with the perturbations in the cold dark matter, they are 
not as significant as can be observed from the evolution of log(^) and log(^) 
vs. log(?7) in figure 3.4. In this case perturbations in a are even less significant 
than in the radiation case. We note that the fact that the linear regime is a very 
good approximation at early times, since as can be seen from figure 3.4, ^ tracks 
oc i^/^ at early times. From the detail of figure 3.4 we see that variations of the 
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Figure 3.3: The evolution of a{ri) and a{ri) for dust- dominated universes with k — Q 
(dashed) and k — 1 (solid). 





Figure 3.4: The evolution of log ^ (solid) and log (dashed) vs log{r]) for a 
dust- dominated universe. The dotted lines correspond to evolution in the linear 
perturbation solution. 
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cold dark matter will start to occur before than variations of a, and ^ is at least 
three orders of magnitude bigger than This is an indication that perturbations 
in the matter fields will drive perturbations in a. 

3.6 Summary and Comments 

The results presented in this chapter, are a quite good approximations when we con- 
sider large scales, but are expected to break down when extended to small scales, 
where non-linear effects come into play and local deviations from isotropy and homo- 
geneity of the matter content are significant. We note that the background solutions 
for ■0, about which we have linearised the perturbations of the Einstein equations, 
are solutions which describe the time evolution of t/^ on a 'standard' FRW back- 
ground. Our neglect of the back-reaction of the ip perturbations on the background 
expansion dynamics is a good approximation up to logarithmic corrections. 

When we examined the non-linear evolution of spherical inhomogeneities by 
means of a comparative numerical study of fiat and closed Friedmann models we 
found that perturbations in a remain almost negligible with respect to perturbations 
in the fluid that dominates the energy density of the universe at the same epoch. 
Comparing the flat with the closed solution for a, we concluded that in both the 
cases of a radiation or dust-dominated epoch, a will 'feel' the change of behaviour 
in the scale factor, and the perturbations ^ will grow in time. The early-time be- 
haviour of the non-linear solutions conflrms the linear behaviour found previously. 
In particular ^ changes in proportion to ^ at early times. 

We have provided a detailed analysis of the behaviour of inhomogeneous per- 
turbations in a and its time variation on large scales under the assumption that 
the defining constant of the BSBM theory, C, is constant and negative in sign. In 
reality this assumption will break down on small scales. The negativity of the ef- 
fective value of ( requires that the cold dark matter is dominated by the magnetic 
rather than the electric field energy (see also the discussion of reference [107] and by 
Bekenstein [83]). However, on sufficiently small scales the dark matter will become 
dominated by baryons and the sign of the effective C will have to change sign. Over- 
all, there will also be a gradient in the value of |C| refiecting the scale dependence 
of the relative contribution of dark matter to the total density of the universe. We 
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have not included these effects in the present analysis. They would need to be 
included in any detailed analysis of the small-scale behaviour of inhomogeneities in 
a. This is an important challenge for future work because it would enable the quasar 
data on varying a to be compared directly with the limits from the Oklo natural 
reactor [8, 22] and Rhenium- Osmium abundances in meteorites [25, 21]. At present 
the relation between the cosmological and geonuclear evidences is unclear because 
the latter are derived from physical processes occurring within the cosmologically 
non-evolving solar system environment. 

Variations of a also affect the cosmic microwave background radiation spectrum 
and anisotropy in different ways, but the effects must be disentangled from allowed 
changes in other cosmological parameters that can contribute similar effects. These 
changes in the microwave background with a left as a free constant parameter were 
analysed in [46] using the new WMAP data. They are far less sensitive that the 
many-multiplet analyses of quasars at z = 0.5 — 3.5 [35, 5, 31, 32], although they 
derive from higher redshifts, z < 1100. These studies can accommodate constant 
and varying a but up to a level that would be too large to be consistent with the 
quasar data and the slow time-evolution of the theory with time- varying a described 
in this thesis. 



CHAPTER 4 



Spherical Collapse Model in the Presence of 

Varying Alpha 



...e neste horizonte do instante que se espelha a liberdade, 
um rosto sem fundo sob um silencio mudo: o eco da temporalidade. 

- Ernesto Mota - 



4.1 Introduction 

In the last chapter, using the hnear theory of the cosmological perturbations, it was 
shown that small inhomogeneities in the matter fields will induce small perturbations 
in the fine structure constant. Those results just confirmed our initial intuition 
that one cannot neglect the effects of gravity in the cosmological evolution of the 
fine structure constant. The reason is simple: any varying-a theory implies the 
existence of a field responsible for variations of the fine structure constant. This 
field is necessarily coupled to the matter fields, at least, to the ones which interact 
electromagnetically. Variations of a, allowed by the conservation of energy and 
momentum will then depend on the expansion of the universe and the evolution 
of the electromagnetically coupled matter. Any inhomogeneities of the later will 
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then affect the evolution of a. In particular, spatial variations of the matter fields 
will then produce spatial variations of the fine structure constant. In reality, we 
will see in this chapter that not only the spatial homogeneity of a is affected by 
the growth of inhomogeneities in matter, but also its time evolution will change 
according to it. This effect is model independent since it is just a product of the 
coupling between the field responsible for the variations of a and matter. The 
inclusion of electroweak or grand unification theories will create even more general 
couplings, since the field, responsible for the a-variations, may also be couple to all 
or some of the other matter fields, besides the ones that interact electromagnetically. 
Due to this, additional consequences will then result [111, 61]. 

It is useful to recall once again the present observational and experimental con- 
straints that a successful varying-a model has to satisfy. We saw in chapter 1, that 
these constraints can be divided into two main groups: local and astro-cosmological. 
The local constraints derive from experiments in our local bound gravitational sys- 
tem: the Oklo natural reactor {z = 0.14), where |^| < 10^^ [8]; the intra solar- 
system decay rate ^^^Re ^^^Os, {z = 0.45), where |^| < 10"^ [57]; and the 
stability of terrestrial atomic clocks {z — 0), where |^| < 4.2 x 10~^^ yr~^[15]. 
Other limits arise from weak equivalence experiments [42] but the limits they pro- 
vide are more model dependent. The astro-cosmological constraints are: the quasar 
absorption spectra z — 0.5 — 3.5, where |^| ^ 10~^; the cosmic microwave back- 
ground radiation {z — 10^), where |^| < 10~^ [46] and Big Bang nucleosynthesis 
{z = 10« - 10^°), where |^| < 2 x 10"^ [46]. 

There is a potential discrepancy between local and astro-cosmological constraints; 
in particular, between the constraint of Aa/a < 10~^ at redshift z — 0.45, com- 
ing from the /3-decay in meteoritic samples [57], and the explicit variation in a of 
Aa/a 10~^ at 2; = 0.5, coming from the low-redshift end of the quasar absorption 
spectra [5] . A successful theory of varying a needs to explain this difference. This is 
a challenge. Models which use a very light scalar field, to drive variations in a, need 
an extreme fine tuning in order to satisfy the phenomenological constraints coming 
from geochemical data (Oklo, /3-decay), the present equivalent principle tests, and 
the quasar absorption spectra, simultaneously [22]. 

We will see in this chapter, that just by looking at the evolution of a in a uni- 
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verse where non-linear structures are formed, the discrepancy between geochemical 
constraints and the quasar data becomes "naturally" justified. 

4.2 The Spherical Collapse Model 

In order to study the behaviour of the fine structure constant during the formation 
of non-hnear overdensities we will use the spherical infall model [128, 52, 157]. This 
will give us an idea on how the evolution of a in the overdensity "breaks away" from 
the corresponding one in the background. 

An overdense sphere is a very useful non-linear model, as it behaves in exactly 
the same way as a closed sub-universe. The density perturbations need not to be a 
uniform sphere: any spherically symmetric perturbation will evolve at a given radius 
in the same way as a uniform sphere containing the same amount of mass [128] . In 
what follows, therefore, density refers to mean density inside a given sphere. 

Let us consider a spherical perturbation with constant density inside it which, 
at an initial time, has an amplitude 5j > and \5i\ ^ 1. For the case of a flat 
matter-dominated FRW universe it is possible to find an exact solution for the 
evolution of the overdensity, although that is not the case for more complex cases, 
like ours, and a numerical simulation is required. Nevertheless the general behaviour 
of the overdensity will be similar. At early times the sphere expands together with 
the background. For a sufficiently large 5i gravity prevents the sphere to expand 
forever. In that case, there are three interesting epochs in the final stages of its 
development. Turnaround: The sphere breaks away from the general expansion 
and reaches a maximum radius. Collapse: If only gravity operates, then the sphere 
will collapse to a singularity where the densities of the matter fields would go to 
infinity. In practise the collapse will never happen. As overdensities break away from 
the background expansion, Virialisation occurs, and they become gravitationally 
bounded systems through a process of violent relaxation [158, 159]. The underlying 
mechanism of violent relaxation is the chaotic gravitational field of a collapsing, 
non-spherical self-gravitational system. The time and space varying gravitational 
field provides the means for individual, non-interacting particles to change their 
energies and to become well mixed in phase space. After only a few dynamical times, 
(r ~ {Gp)~^^'^) the result is a viriahsed distribution of matter whose phase space 
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distribution is roughly Maxwellian and whose density varies as ~ [128]. Due to 
its spherical symmetry and isothermal solution, such a configuration is referred as an 
isothermal sphere. In addition, baryonic matter can, through dissipative processes, 
lose energy and condense even further into cores of the isothermal spheres that form. 
The process of structure formation is much more complex and we shall not attempt 
to describe it here, referring the reader to [52, 128]. 

With virialisation, a condition of equilibrium is then reached and the system 
becomes stationary: There are no more variations of the radius of the system, or in 
its energy components. The stationary condition reached by the process of virialisa- 
tion is very important, and affects dramatically the variations of the fine structure 
constant, as we will see bellow. 

Once again we will use the BSBM theory in order to study the effects of structure 
formation in the evolution of a. We will assume a A-Cold-Dark-Matter [AC DM) 
model of structure formation. 

The background evolution of the universe and its components is given by the 
equations (1.7.7), (1.7.8) and equation (1.7.9). This equations will describe the 
evolution of the background scale factor, a, the scalar field, tp, responsible for the 
variations of a and the energy density of non-relativistic matter p^, respectively. 

Just to clarify to notation from now on, we will label some variables which "live 
in" the background with a b. Quantities which "live" inside the overdensities will 
be labelled by c. For instance, the value of a inside a cluster will be referred to etc. 
Quantities which represent the epoch at or after virialisation, will be labelled by a v. 
For instance, Aay/a^, represents the "time density contrast" at virialisation. This 
quantity represents a comparison between the value of the fine structure constant 
in a virialised clusters or in the background, at a given redshift, and its value on 
Earth, ag. 

4.2.1 Dynamics of the overdensities 

Variations of a inside overdense spherical regions can be studied using the spherical 
infall model. The evolution of a spherical overdense patch of scale radius R{t) is 
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given by the Priedmann acceleration equation: 

I = {Pcam (1 + Id e-^'^=) + 4p^. - 2pa) (4.2.1) 

where Pcdm is the total density of cold dark matter and baryons in the cluster, ilj^. is 
the homogeneous field inside the cluster, pA is the cosmological constant density and 
p^^ = We have used the equations of state p^^ — p^^, Pcdm — and pa = ~P\-i 

where p represents the pressure of the fluid. 

In the cluster, the evolution of '0c and pcdm are given by 

4 + 3^4 = --e-^^%Pcdm (4.2.2) 
R 

Pcdm = —^^Pcdm (4.2.3) 

According to the Virial theorem, a condition of equilibrium will be reached when 

1 dU , , 

where, T is the average total kinetic energy and U is the average total potential 
energy in the sphere. Note that we obtain the usual T — condition, for any 
potential with a powerlaw form {U oc R^), which will be our case. 

The potential energy for a given component x can be calculated from its general 
expression in a spherical region [160]: 

C/^ = 27r / ptot^^dr (4.2.5) 



</.,(r) = -2T:Gp:,[R'--j (4.2.6) 

where ptot is the total energy density inside the sphere, (f)x is the gravitational field 
potential due to the p^ density component. 
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In the case of a AC DM model the potential energies inside the cluster would be: 

Ug = -IcM^R-^ (4.2.7) 
5 

Ua = -^tiGpaMR^ (4.2.8) 
5 

U^^ = —GMM^^R-^ (4.2.9) 

where Ug is the potential energy associated with the uniform spherical overdensity, 
Ua is the potential associated with A, and U^^ is the potential associated with -^c- 
M — Mcdm + M/Jc is the cluster mass, with 

M,dm = —pcdmil + \C\e-^'^')R^ (4.2.10) 
= yP^c^' (4.2.11) 



The virial theorem would be satisfied when: 

Tyir ^—Ug + Ua- 2U^^ (4.2.12) 

where Tyir = jMv^^j. is the total kinetic energy at virialisation and v'^^^ is the mean- 
square velocity of masses in the cluster. 

Using the virial theorem (4.2.12) and energy conservation at the turn around 
and when the cluster virialises, we obtain an equilibrium condition in term of the 
potential energies 

^UGiZy) + 2C/a(^„) - U^^iz^) = UGizta) + UAizta) + U^^a) (4.2.13) 

where is the rcdshift of virialisation and Zta is the redshift at the turn-around 
of the over-density at its maximum radius, when R = Rmax and i? = 0. In the 
case where C = and ^ = ^/^ = we just have the usual virialisation condition for 
AC DM models [161, 162]. 

The inclusion of a varying a at a level consistent with observation [35, 5, 30, 31, 
32, 33], {a/ao ~ 10~^H), does not affect the overall expansion of the universe up to 
logarithmic corrections, if we are far from the initial singularity [108]. The reason 
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for that is the neghgible contribution of to the energy content of the universe, 
see figure 1.2. In a similar way, the dynamical collapse of the overdense regions 
is also not affected: The energy density associated with is always a neghgible 
contribution to the energy content of cluster. Even the kinetic energy, p^^ — ^ip'^, 
which starts to grow in the final stages of the cluster evolution will be negligible, 
see Figure 4.1. 




log a 

Figure 4.1: Evolution of log pdm (solid lines) , logp^^ (dash- dotted line), logpA 
(dotted line) and log pa (dashed line) inside an overdensity, as a function oflog{a), 
in a ACDM model. Where pa = p^c + \C\Pcdme~'^'^'' ■ 



4.3 Evolution of a During the Cluster Formation 

The behaviour of the fine structure constant during the evolution of a cluster can 
now be obtained numerically evolving the background equations (1.7.7)-(1.7.9) and 
the cluster equations (4.2.1)-(4.2.3) until the virialisation condition holds, equation 
(4.2.13). 
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4.3.1 Setting the initial conditions 

In order to satisfy the constraints imposed by the observations, we need to set up 
initial conditions for the evolution. Since the Earth is now a virialised overdense 
region, the initial condition for is chosen so as to obtain our measured laboratory 
value of a at virialisation, ac{zv) = = ao- Since the redshift at which our 
cluster has virialised is uncertain, we will choose a representative example where 
virialisation occurs over the range < Zy < 10. 

After we have set for Earth, another constraint we need to satisfy is given by 
the quasar observations [35]. This means that when comparing the value of the fine 
structure constant on Earth, at its virialisation, = ao, with the value of the fine 
structure constant of another region at some given redshift in the range accessed by 
the quasar spectra, 3.5 > z > 0.5, we need to obtain Aa/a = {a{z) — q;„)/q;„ ~ 
—5.4 X 10~^. This raises the question as to the location of the clouds where the 
quasar absorption lines are formed: are they in a region which should be consider 
as part of the background or an overdensity with somewhat lower contrast than 
exists in our Galaxy? Unfortunately, this question cannot be answered because we 
do not know the density of the clouds, only the column density. Nevertheless, these 
clouds are much less dense than the solar system. Because of this, it is a good 
approximation to assume that the clouds possess the background density. 

Thus, the initial conditions for ip are chosen so as to obtain our measured labo- 
ratory value of a at virialisation a^Zy) — ckq and to match the latest observations 
[35] for background regions at 3.5 > \z — z^] > 0.5. 

4.3.2 Differences between the background and the overden- 
sities 

Figure 4.2 show us the evolution of the fine structure constant inside an evolving 
overdensity and in the background. It is clear that there will be differences between 
the value of the fine structure constant inside a cluster, etc, and in the background, 
ai). It is also shown in that figure, that the value of a will also dependent on the 
redshift at which an overdensity will virialise. 

In figure 4.3, we have evolved four clusters, which have different initial conditions 
in order to satisfy the constraints described in the previous section. This is just an 
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Figure 4.2: Evolution of a in the background (dashed line) and inside clusters (solid 
lines) as a function o/log(l + z). Initial conditions were set to match observations of a 
variation of ref. [5]. Two clusters virialise at different redshifts, one of them in order to 
have c{zv = 1) = oeo- Vertical lines represent the moment of turn-around. 
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Figure 4.3: Evolution of Aa/a in the background (dashed lines) and inside clusters (solid 
lines) as a function o/log(l + z). Initial conditions were set to match observations of a 
variation of ref. [5]. Four clusters that virialised at different redshifts. All clusters were 
started so as to have ac{zv) = ao- 
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example, since in reality, the initial condition for needs to be fixed only once, 
for our Galaxy. Hence, a in other clusters will have a lower or higher value (with 
respect to ckq) depending on their Zy] as can be seen in figure 4.2. The example given 
was just to show that independently of the redshift at which Earth virialises, the 
order of magnitude between virialised regions and the background is always of order 
10~^. Of course this is not a coincidence, but a consequence of setting the initial 
conditions in order to satisfy the latest quasar observations [35] . This justifies why 
our study, although is a qualitative analysis, also can give us a good approximation 
of the order of magnitude expected for variations of the fine structure constant both 
inside overdensities and in the background. The choice of having assumed that the 
clouds, where the quasar observations were made, correspond to the background, is 
then a not very stringent condition. 

4.3.3 Constancy of a after virialisation 

After virialisation occurs, the cluster radius, i?, becomes constant; time and space 
variations of a are suppressed, and a becomes constant. If there were any variations 
of a after virialisation, the energy components and radius of the cluster would need 
to vary in order to conserve energy and momentum. This would inconsistent with 
virialisation. This phenomenon is not included in figure 4.2, since we did follow the 
evolution to virialisation with a many-body simulation which would need to include 
the fluid equations that describe the pressure inside the cluster, see for instance 
[163, 164]. In our simulation, the virial condition is a 'stop condition' and so we just 
observe the typical behaviour of the cluster's collapse as i? — > 0, see figure 4.4, or 
Pcdm +00, see figure 4.1. 

In addition to the stationarity condition that occurs when the cluster virialises, 
it can be seen from figure 4.3 that in all cases the variation of a since the beginning 
of the cluster formation is of order 10^^, and numerical simulations give a /a ~ 
10~^^s~^, we will see this in the next chapter. If variations of a are so small for such 
a wide range of virialisation rcdshifts, we can assume that the difference between the 
value of «c at and at 2; = will be negligible. Therefore it is a good approximation 
to assume that the time-evolution of both etc and of the cluster will cease after 
virialisation. Although this is not necessarily true ( for instance the cluster could 
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Figure 4.4: Evolution of the radius of a cluster, R, as a function o/log(l + z). 

keep accreting mass), it is a good approximation in respect of the evolution of a, 
especially for clusters which have virialised at lower redshifts. 

The fact that local a values 'freeze in' at virialisation, means we would observe 
no time or spatial variations of a on Earth, or elsewhere in our Galaxy, even though 
time-variations in a might still be occurring on cxtragalactic scales. For a cluster, 
the value of a today will be the value of a at the virialisation time of the clus- 
ter. We should observe significant differences in a only when comparing clusters 
which virialised at quite different redshifts. Differences will arise within the same 
bound system only if it has not reached viral equilibrium. Hence, variations of a 
using geochemical methods could easily give a value that is 10 — 100 times smaller 
than is inferred from quasar spectra. For instance, in spite of the fact that we 
cannot simulate the evolution of etc after virialisation, it is useful to compare the 
value of Aa/a in clusters that have have virialised at a low redshift and in the 
background, to have an order of magnitude of the variations in the fine structure 
constant. Numerical simulations, normalised to have adzy — 0) — ao and to satisfy 
the quasar observations, give Aa/a{z — 2) — —6 x 10~^ for the background and 
Aa/a{zv — 0.13) = —1.5 x 10~^ in a cluster which has virialised at Zy — 0.13. 
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4.4 Measuring a in Virialised Regions 

If we could measure the fine structure constant inside viriafised overdensities and 
the corresponding value in the background, at the time of their virialisation, then 
figure 4.5 would give us the evolution of or its time shift, Aay/ay, as a function 
of z^. 




>,7.297335 - 



to 7.297330 - 
O 

X 7.297325 - 
7.297320 - 

a 

"7.297315 - 



7.297310 - 




log ( 



Figure 4.5: Plot of a and Aa/a as a function o/ log(l + z), at virialisation. Clusters 
(solid line), background (dashed line). 
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Those figures show us that differences in a between the background and the 
overdensities increase as 2;„ — > 0. This is due to the earher freezing of the value of a 
at viriahsation, and to our assumption that we hve in a ACDM universe. At lower 
redshifts, specially after A starts to dominate, variations of a in the background are 
turned off by the accelerated expansion. However, the value of a in the collapsing 
cluster will keep growing until viriahsation occurs. At higher redshifts, z ^ 1, both 
ab and etc evolve in expanding environments: their increase is logarithmic in time 
before A starts to dominate, so the difference between them will be much smaller. 

These results mean that if we were able to compare the value of a, in overdensities 
which have virialised at a similar redshift as ours, with regions that have virialised 
much before or much after, we would notice that there were indeed differences in 
the measured Aa/a. 

4.4.1 Dependence on the matter density of the clusters 

Variations of etc with respect to ab can be tracked using the "spatial density con- 
trast" , 

Sa _ ac — ab 

a ab ' 

and computed at viriahsation or while the cluster is still evolving. 

A plot of the "spatial density contrast" in a with respect to the matter density 
inhomogeneity, Sp/p, at viriahsation is shown in figure. 4.6. Note that 5a/ a grows 
in proportional to the density contrast of the matter inhomogeneities. 

In a ACDM model, the density contrast, Ac = Pcdmi^v)/ Pbi^v) increases as 
the redshift decreases [161]. For high redshifts, the density contrast at viriahsation 
becomes the asymptotically constant in standard (A = 0) CDM, Ac ^ 178, if pb is 
measured at collapse or Ac ~ 148, if pb is measured at viriahsation. This is another 
reason why at lower redshifts the difference between a^ and ab increases. 

Trends of variation of a can be then predicted from the value of the matter 
density contrast of the regions observed. Useful fitting formulae for the dependence 
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Figure 4.6: Variation of 5a/a with 5p/ p at the cluster virialisation redshift. The evolution 
of a inside the clusters was normalised to satisfy the latest time-variation observational 
results, and to have adzy) = ao for Zy= 1. 

of a variation on Sp/p and the scale factor a = (1 + z)^^ are: 

Sex 1 3 

_ = (5.56 - 0.7(75 + 0.078(7 + 0.00352(7^) x 10"^ 

a 

— = (5.37 + 0.373^5 - 0.27^3 + 0.007^5) x 10"^ (4.4.1) 
a 

— = -(2.47-1.81(75 +0.59(7- 0.094(75) X 10"^ 
a 

where 9 = ^/Ac„, a = a /ay. Ac„ and are "input" parameters defined by our 
choice of the redshift of when a{zy) — ao . These formulae are valid for clusters 
which virialised at < 2;„ < 15. For the case where a{zy — 1) — ao these formulae 
are accurate to better than 1%. 

4.5 Comments and Discussion 

Generally, past studies of spatially homogeneous cosmologies have matched the value 
of aj, with the terrestrial value of a measured today. However, it is clear that 
the value of the fine structure constant on Earth, and most probably in our local 
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cluster, will differ from that in the background universe. This feature has been 
ignored when comparing observations of a variations from quasar absorption spectra 
[35, 5, 31] with local bounds derived by direct measurement [165, 166] or from Oklo 
[20, 8, 22] and long-hved /3-decayers [57]. A similar unwarranted assumption is 
generally made when using solar system tests of general relativity to bound possible 
time variations of G in Brans-Dicke theory [167]: there is no reason why G/G should 
be the same in the solar system and on cosmological scales. Since any varying- 
constant model require the existence of a scalar field coupled to the matter fields, 
our considerations apply to all other theories besides BSBM and to variations of 
other "constants" [106, 7]. The consideration of the inhomogeneous evolution of the 
universe is therefore essential for a correct comparison of extra-galactic and solar 
system limits on, and observations of, possible time variation of the fine structure 
constant and other "constants" . 

There are many different constraints, a theoretical model of varying-a, needs 
to satisfy, in order to be in agreement with the observations and experiments. In 
particular, it needs to satisfy the geochemical constraints at, z = 0.45, of l^aja < 
10~^ and the quasar data, at z — 0.5, of Aa/a ~ 10"^. 

In order to explain, why any experiment made on earth [20, 8], or in our local 
system [167, 21], present constraints on variations of a which are, 10 to 100 times 
lower than the quasar spectra observations [35], we need to take into account the 
mechanism of structure formation in the universe. We have then considered that 
today's value of a is measured on earth (and on our solar system), and that these 
are virialised overdensities. We have matched the quasar absorption spectra obser- 
vations [35], comparing the today's value of a on a virialised region (our galaxy) and 
the background (which represented the clouds where the quasar data was obtained) . 
Using these assumption, to set up our initial conditions, we have used the spherical 
collapse model to investigate the effects of non-linear structure formation on the 
cosmological evolution of the fine structure constant. Evolving a tiny overdensity 
from the radiation era up to today, we were able to observe the differences between 
the evolution of a inside the overdensity and the background. When the density of 
these structures becomes big enough for them to break away from the background 
expansion, they start to collapse until its virialisation. When the overdensity viri- 
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alises, variations of a are highly suppressed, and it becomes constant inside that 
region, otherwise the system would not be able to virialise. 

The 'freeze in' effect of the fine structure constant after virialisation is of major 
importance. Any experiment or observation made, in order to study variations 
of a, on Earth or in our local system , would measure no variation at all. Or if 
they do, those would be much smaller than the results obtained from the quasar 
absorption spectra. So, variations of a, would naturally satisfy any experiments on 
the violation of the principle of equivalence [56, 57, 58] , in our local overdensity, and 
any geochemical observation. At the same time we would be able to explain the 
quasars absorption spectra observations [35, 5, 30, 31]. 

The dependence of a on the matter-field perturbations is much less important 
when one is studying effects on the early universe, for example on the CMBR or Big 
Bang nucleosynthesis [46]. In the linear regime of the cosmological perturbations, 
small perturbations in a will decay or become constant in the radiation era [2]. Also, 
in the non-linear regime of structure formation, the evolution of the fine structure 
constant in the background and inside the overdensities will not differ until the 
overdensities turnaround. Turnaround occurs at a far lower redshift than that of 
last scattering, z = 1100. At these high redshifts, it was found that Aa/a < 10~^. In 
the background, the fine structure constant, will be a constant during the radiation 
epoch [107]. A growth in value of a will happen only in the matter dominated era 
[1]. Hence, the early-universe constraint, coming from the CMBR and Big Bang 
Nucleosynthesis, Aa/a < 10^^, are easily satisfied. 



CHAPTER 5 



The Fine Structure Constant Dependence on the 
Dark-Energy Equation of State and on the 

Coupling to Matter 



O todo 6 tudo, e o mini que sou e o resto, 
ndo hd resto, hd uma angustia isso sim 
de saber que o todo que sou ndo o set dizer... 

mas o sei que sou todo assim... 

- Ernesto Mota - 



5.1 Introduction 

In the last chapter, we saw that the cosmological evolution of the fine structure 
constant depends on the evolution of non-linear structures in the universe. There 
are two main features which affect the evolution of o; in a universe where large scale 
structures are formed. The most obvious, is the coupling between the scalar field, 
which drives variations in a, and the matter fields. The second is the dependence 
of non-linear models of structure formation on the equation of state of the universe, 
and in particular that of the dark energy one [161, 162]. 
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Looking at the right-hand-side of equation (1.7.8), we see that, C,/uj represents 
the couphng between the scalar field ^ an the matter density. It is then natural to 
think that the rate of change of the fine structure constant, both in the background 
and inside overdensities, will be a function of the absolute value of C,/^- Smaller 
values of |C/<^| would lead to a smaller dependence of the variation of a on the 
matter inhomogeneities. 

Besides the dependence on the coupling C/o;, the fine structure constant also 
depends on the equation of state of the universe. For instance, several solutions 
were found for cct, for the different epochs the universe went through. In the case of 
the overdensities, the dependence of etc on the equation of state of the universe is 
related to the dependence of the density contrast of the overdensities on the structure 
formation model. 

In resume, the evolution of the fine structure constant in the background, and 
inside clusters then depends mainly on the dominant equation of state of the universe 
and the value and sign of the coupling constant ^. We will then investigate the 
inhomogeneous evolution of the fine structure constant after the development of 
non-hnear cosmic structures, that is, after they have virialised. This allow us to 
predict variations of the fine structure constant among virialised systems and to 
investigate the dependence of a on the equation of state of the dark energy. In the 
last section, we show how spatial variations in a may occur due to possible spatial 
variations in the coupling of a to the matter fields. 

5.2 BSBM model and Dark Energy Structure For- 
mation Models 

In order to study the dependence of the evolution of a on the equation of state of 
dark energy, we need to generalise the equations used in the previous chapter both 
for the background and in the overdensities. 

The background universe will be described by a fiat, homogeneous and isotropic 
Priedmann metric with expansion scale factor a{t). The universe contains pressure- 
free matter, of density pm oc and a dark energy fiuid with a constant equation 
of state, ty^, and an energy-density oc a"^^^"*'"'*). In the case where dark energy 
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is the cosmological constant A, then = p\ = A/{87iG) and — —1. Hence the 
Priedmann equation (1.7.7) becomes: 

H'^^ {pm (1 + Id e-'^) + + p^) (5.2.1) 

The evolution of a spherical ovcrdcnse patch of scale radius R{t) is the given by 
the Friedmann acceleration equation (4.2.1), which becomes: 

I = {Pcdm (1 + Id e-'^O + 4p^. + (1 + 3«;^Jp^J (5.2.2) 

where pcdm is the density of cold dark matter in the cluster, p<^^ is the energy density 
of the dark energy inside the cluster and p^^ = l'^^, where ■^c represents the scalar 
field inside the overdensity. We have also used the equations of state p^^ — p^^, 
Pcdm = and p^^ = w^^p^^. 

In the cluster, the evolution of ■0c, Pcdm and p^^ are given by 



R 

i'c + 3— v^c = 


--e-^Hpcdm 

UJ 


(5.2.3) 


Pcdm — 


-.^ 

'-> jyPcdm 


(5.2.4) 


P^c = 


-3^(l+W^c)P<Ac 


(5.2.5) 



Once again wc will evolve our overdensity from a redshift well inside the radiation 
era until its virialisation occurs. In order to know when the cluster becomes virialised 
we will use the Virial theorem, equation (4.2.4). Again it will be useful to write the 
condition for virialisation to occur, using the potential energies associated to the 
different components of the overdensity. 

Since we are now considering a dark energy fluid with a general, but constant, 
equation of state w^^, we need to generalise the potential energy associated to the 
cosmological constant, U\. The generalisation of equation (4.2.9), in order to incor- 
porate a fluid with a general equation of state is trivial. Using equations (4.2.6) and 
(4.2.9), we obtain the potential energy associated to the dark energy fluid 

U^^ = ~GM^MR'<'+^^^) (5.2.6) 
o 
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where M^^ is the mass associated to that fluid, given by 

Using now U^^ in equation (4.2.4), instead of [/a we obtain a condition for viri- 
ahsation, similar to equation (4.2.13), in terms of the known variables. 

The behaviour of the flne structure constant during the evolution of a cluster can 
now be obtained numerically evolving the background equations (1.7.8) and (5.2.1) 
and the cluster equations (5.2.2)-(5.2.5) until the viriahsation holds. 

Once again the initial conditions for are chosen so as to obtain our measured 
laboratory value of a at viriahsation ac{zy) = ay — ag and to match the latest 
observations [35]: Aa/a fa —5.4 x 10""^ for regions at 3.5 > l-z — -z^,! > 0.5, where Zy 
will be the redshift at which our inhomogeneity virialised. 

5.3 Dependence on the Dark-Energy Equation of 

Sate 

Non-linear models of structure formation will present different features depending 
on the equation of state of the universe [161, 162]. The main difference is the way 
the parameter Ac = Pcdmi^c) / Pbi^v) evolves with the redshift. The evolution of Ac 
depends on the equation of state of the dark-energy component which dominates 
the expansion dynamics. For instance, in a ACDM model, the density contrast, 
Ac increases as the redshift decreases. At high redshifts, the density contrast at 
viriahsation becomes asymptotically constant in standard (A = 0) CDM, with 
Ac ~ 178 at collapse or Ac ~ 148 at viriahsation. This behaviour is common 
to other dark-energy models of structure formation {wCDM), where the major 
difference is in the magnitude and the rate of change of Ac at low redshifts. 

The local value of the fine structure constant will be a function of the redshift and 
will be dependent on the density of the region of the universe we measure it, accord- 
ing to whether it is in the background or an ovcrdcnsity [3]. The density contrast of 
the virialised clusters depends on the dark-energy equation of state parameter, w,^ . 
Hence, the evolution of a will be dependent on w</, as well. 
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What difference we would expect to see in a if we compared two bound systems, 
like two clusters of galaxies? And how does the difference depend on the cosmological 
model of structure formation? These questions can be answered by looking at the 
time and spatial variations of a at the time of virialisation. The space variations 
will be tracked using a "spatial" density contrast, 

— = — (5.3.1) 

a ah 

which is computed at virialisation (where ac{z = z^) = otv). 

Since the main dependence of a is on the density of the clusters and the redshift 
of virialisation, we will only study dark-energy models where W(j, is a constant. Any 
effect contributed by a time- varying equation of state should be negligible, since the 
important feature is the average equation of state of the universe. This may not be 
the case for the models where the scalar field responsible for the variations in a is 
coupled to dark energy [168, 169, 170]. 

In order to have a qualitative behaviour of the evolution of the fine structure con- 
stant, at the virialisation of an overdensity, we will then compare the standard Cold 
Dark Matter model, SCDM, to the dark-energy Cold Dark Matter, wCDM, mod- 
els. In particular, we will examine the representative cases of W(j) — —1, —0.8, —0.6. 
All models will be normalised to have — ao at z — and to satisfy the quasar ob- 
servations [35], as discussed above. This normalisation, although unrealistic (Earth 
did not virialisc today), give us some indication of the dependence of the time and 
spatial evolution in a on the different models. In reality, this approximation will not 
affect the order of magnitude of the spatial and time variations in a for the cases of 
virialisation at low redshift [3]. 

5.3.1 Time-shifts and the evolution of a 

The final value of a inside virialised overdensities and its evolution in the background 
is shown in figure 5.1. Prom this plot, a feature common to all the models stands out: 
the fine structure constant in the background regions has a lower value than inside 
the viriahsed overdensities. Also, its eventual local value depends on the redshift at 
which the overdensity virialises. 
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Figure 5.1: Value of a in the background and inside clusters as a function ci/log(l + 
Zy), the epoch of virialisation. Solid line corresponds to the standard (A = CDM, 
dashed line is the ACDM, dashed-dotted corresponds to a wCDM model with = 
—0.8, dotted-line corresponds to a wCDM model with = —0.6. In all the models, 
the initial condition were set in order to have a^z — 0) — ag and to satisfy Aa/a ~ 
-5.4 X 10-6 at 3.5 > U - zj > 0.5. 
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As expected, the equation of state of the dark energy affects the evolution of 
a, both in the overdensities and the background. A major difference arises if we 
compare the SCDM and wCDM models. In a SCDM model, the fine structure 
constant is always a growing function, both in the background and inside the over- 
densities, and the growth rate is almost constant. In a wCDM model, the evolution 
of a will depend strongly on whether one is inside a cluster or in the background. 
In the wCDM background, ctb becomes constant (independent of the redshift) as 
the universe enters the phase of accelerated expansion. Inside the clusters, will 
always grow and its value now will depend on the redshift at which virialisation oc- 
curred. The cumulative effect of this growth increases as we consider overdensities 
which virialise at increasingly lower redshift. 

These differences arise due to the dependence of the fine structure constant on the 
equation of state of the universe and the density of the regions we are measuring it. 
In a SCDM model, we will always live in a dust-dominated era. The fine structure 
constant will then be an ever-increasing logarithmic function of time, a oc ln(t) 
[107, 1]. The growth rates of and Uc will be constant, since Ac is independent of 
the redshift in a SCDM model. In a wCDM model, dark energy plays an important 
role at low redshifts. As wc reach low redshifts, where dark energy dominates the 
universal expansion, ah becomes a constant [107, 1], but Ac continues to increase, 
as will ttc [3]. The growth of a^, becomes steeper as we go from a dark-energy fluid 
with Wfj) = —0.6 to the A-likc case of w^f, = —1. The intermediate situation is where 

= —0.8, due to the dependence of Ac on w^. 

Similar conclusions can be drawn with respect to the "time shift" of the fine 
structure constant, (Aa/a), at virialisation, see figure 5.2 
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Figure 5.2: Variation of Aa/a in the background and inside clusters as a function 
of the epoch of virialisation, log(l + Zy). The solid line corresponds to the standard 
(A = 0) CDM, the dashed line is the AC DM, the dashed-dotted line corresponds 
to a wCDM model with = —0.8, the dotted line corresponds to a wCDM model 
with — —0.6. In all models, the initial conditions were set in order to have 
ac{z — 0) — ao and to satisfy Aa/a ~ —5.4 x 10~^ at 3.5 > \z — Zy\ > 0.5. 
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Figure 5.3: Variation of 5a/ a as a function o/log(l + z^). Solid line corresponds 
to the standard (h. = Q) CDM, dashed line is the ACDM, dashed-dotted line cor- 
responds to a wCDM model with = —0.8, dotted line corresponds to a wCDM 
model with = —0.6. In all the models, the initial conditions were set in order to 
have ac{z — 0) — and to satisfy Aa/a ~ —5.4 x 10"^ at 3.5 '>\z — Zy\ > 0.5. 



5.3.2 Spatial variations in a 

Spatial variations in a will be dramatically different when comparing the standard 
CDM model with the wCDM models. In the SCDM model, the difference between 
the fine structure constant in a virialised cluster (a^,) and in the background {ai,) 
will always be the same, 5a/ a ^ 5.2 x 10~^, independently of the redshift at which 
we measure it, figure 5.3. Again, this is because, in a SCDM model, Ac is always 
a constant independent of the redshift at which virialisation occurs. The constancy 
of 5a /a is a signature of the SCDM structure-formation model, and it may even 
provide a means to rule out the SCDM model completely if, when comparing the 
value oi5a/a in two different clusters, we do not find the same value, independently 
of A similar result is found for the case of a dark-energy structure formation 
model at high redshifts where 5a/ a will be constant. This behaviour is expected 
since at high redshift any wCDM model is asymptotically equivalent to standard 
CDM. 

As expected, it is at low redshifts that the difference between wCDM and 
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Figure 5.4: Variation of 5a /a as a function as a function of Ac- The single-point 
corresponds to the standard (A = 0) CDM , the dashed line is the ACDM, the dash- 
dot line corresponds to a wCDM model with = —0.8, dotted-line corresponds to 
a wCDM model with = —0.6. In all the models, the initial conditions were set in 
order to have a^z — 0) — ao and to satisfy Aa/a ~ —5.4 x 10~^ at 3.5 > l-z — -2^1 > 
0.5. 

SCDM emerges. When comparing virialised regions at low redshifts, 5a/ a will 
increase in a wCDM model as we approach z = This is due to an increase of the 
density contrast of the virialised regions, A^, and the approach to a constant value 
of «{,, see figure 5.4. In general, the growth will be steeper for smaller values of w^, 
although there will be parameter degeneracies between the behaviour of different 
models, which ensure that there is no simple relation between 5a/ a and the dark- 
energy equation of state. Note that, independently of the structure formation model 
we use, dic/ac at virialisation is always a decreasing function of time, as shown in 
figure 5.5. 
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Figure 5.5: Variation of logio{a/a) as a function of log {1 + z^). Solid line cor- 
responds to the Standard (A = 0) CDM, dashed-line is the ACDM, dashed-dotted 
corresponds to a wCDM model with w^j) = —0.8, dotted line corresponds to a wCDM 
model with Wtj, — —0.6. In all the models, the initial conditions were set in order to 
have ac{z — 0) — ao and to satisfy Aa/a —5.4 x 10~^ at 3.5 > \z — z^] > 0.5. 
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5.4 Dependence on the Coupling of a to the Mat- 
ter Fields 

The evolution of a in the background and inside clusters depends mainly on the 
dominant equation of state of the universe and the sign of the coupling constant 
C/oj, which is determined by the theory and the dark matter's identity. As was 
shown in chapter 2, will be nearly constant for an accelerated expansion and also 
during the radiation era far from the initial singularity (where the kinetic term, p^, 
can dominate). Slow evolution of a will occur during the dust-dominated epoch, 
where a increases logarithmically in time for C < 0. When C, is negative, a will be 
a slowly growing function of time but a will fall rapidly (even during a curvature- 
dominated era) for C, positive [107]. A similar behaviour is found for the evolution 
of the fine structure constant inside over densities. Thus, we see that a slow change 
in a, cut off by the accelerated expansion at low redshift, that may be required by 
the data, demands that C < in the cosmological background. 

The sign of C, is determined by the physical character of the matter that carries 
electromagnetic coupling. If it is dominated by magnetic energy then < 0, if not 
then C > 0. Baryons will usually have a positive (although Bekenstein has argued 
for negative baryonic C in references [81, 82], but see [22]), in particular ^ 10"*^ 
for neutrons and protons. Dark matter may have negative values of ^, for instance 
superconducting cosmic strings have C ~ — 1- 

In the previous section, we have chosen the sign of C to be negative so a is a 
slowly-growing function in time during the era of dust domination. This was done in 
order to match the latest observations which suggest that a had a smaller value in the 
past [5, 30, 31, 32]. This is a good approximation, since we have been studying the 
cosmological evolution of a during large-scale structure formation, when dark matter 
dominates. However, we know that on sufficiently small scales the dark matter will 
become dominated by a baryonic contribution for which C > 0. The transition in 
the dominant form of total density, from non-baryonic to baryonic as one goes from 
large to small scales requires a significant evolution in the magnitude and sign of 
C/oj. This inhomogeneity will create distinctive behaviours in the evolution of the 
fine structure constant and will not be studied in this thesis. It is clear that a change 
in the sign oiC,/uj will lead to a completely different type of evolution for a, although 
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Figure 5.6: Evolution of a in the background (dashed-line) and inside a cluster 
(Solid line) as a function o/log(l + z). For C/^ = —10^"^, —4 x 10^*^, —7 x 10^'^ and 
— 10~^. The lower curves correspond to a lower value of C,/ui. The initial conditions 
were set in order to have a^z = 0) = ao «^<^ to satisfy Aa/a — —5.4 x 10~^ at 
S > \z — Zy\ > 1, in the case where ( — —7 x 10~^. 

the expected variations in the sign oi C,/uj will occur on scales much smaller than 
those to which we are applying the spherical collapse model here. Hence, we will 
only investigate the effects of changing the absolute value of the couphng, for 
the evolution of the fine structure constant. 

From figure 5.6 it is clear that the rate of changes in and will be functions 
of the absolute value of C,/^- Smaller values of \C/^\ lead to a slower variation of a. 
A similar behaviour is found for the time variations in Aa/a, see figure 5.7. 

The faster variation in a and Aa/a for higher values of \C/^\ is also a common 
feature for 5a/a, see figure 5.8. This is expected. A stronger coupling to the matter 
fields would naturally lead to a stronger dependence on the matter inhomogeneities, 
and in particular on their density contrast, A^. 




5.4 Dependence on the Coupling of a to the Matter Fields 



115 




-9 -8 -7 -6 -5 -4 -3 -2 -1 

log a 



Figure 5.7: Evolution of ^a/a in the background (dashed-line) and inside a cluster 
(Solid line) as a function o/log(l + z). For C,/uj = — 10~^, —4 x 10^'^, —7 x 10^'^ and 
— 10~^. The lower curves correspond to a lower value of (/uj. The initial conditions 
were set in order to have ac{z = 0) = cto and to satisfy Aa/a = —5.4 x 10~® at 
^ > \z — Zy\ > 1, in the case where ( — —7 x 10~^. 
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Figure 5.8: Evolution of 5a /a as a function 5 p/ p. ForC,/uj = — 10""^, — 4x 10""^, — 7x 
10~^ and — 10~^. The lower curves correspond to a lower value of C,/uj. The initial 
conditions were set in order to have ac{z — 0) — ao and to satisfy Aa/a ~ —5.4 x 
10~^ at 3 > \z — Zv\ > 1; in the case where ( — —7 x 10~^. 
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In reality, the dependence of a on the couphng C,/uj has a degeneracy with respect 
to the initial condition chosen for ■0. This is clear from the scale invariance of 
equation (1.7.8) under linear shifts in the value of — > + const and rescaling of 
C,/uj and t. It is always possible, to obtain the same evolution (and rate of change) of 
a and Aa/a for any other value of |C/<^|; see for example Tables 5.1,5.2,5.3 and 5.4, 
where we have tabulated the shifts Aa/a, 5q;/q;, and time change, a/a, obtained 
in the clusters and in the background for various numerical choices oi C,/uj and 
virialisation redshift z^- The observed Oklo and /3-decay constraints on variations 
in a are highlighted in italic and boldface, respectively. 

In the plots of this section, we have chosen the value C,/ijj — —7 x 10~^ to be 
the one which satisfies the current observations, but we could have used any value 
of \C/^\ because of the invariance under rescalings. However, once we set the initial 
condition for a given value of C/o;, any deviation from that value leads to quite 
different future variations in a. This can occur if there are regions of the universe 
where the dominant matter has a different nature, and do possesses a different value 
(and even sign) of (/u to that in our solar system. The evolution of a in those 
regions may then be different from the one that led to the value of a^z = z^) = aQ 
on Earth. 

It is important to note that, due to the degeneracy between the initial condition 
and the coupling to the matter fields, there is no way to avoid evolving differences 
in a variations between the overdensitics and the background. The difference will 
be of the same order of magnitude as the effects indicated by the recent quasar 
absorption-line data. This is not a coincidence and it is related to the fact that 
we have normalised a{z = z^) = ao on Earth and Aa/a = 5.4 x 10~^ at 3.5 > 
1^ — > 0.5. So, any varying-a model that uses these normalisations will create 
a similar difference in a evolution between the overdensitics and the background, 
independently of the coupling (/cu. 

We might ask: how model-independent are these results? In this connection it 
is interesting to note that even with a zero coupling to the matter fields (which is 
unrealistic), C = 0, there is no way to avoid difference arising between the evolution 
of a in the background and in the cluster overdensitics. While the background 
expansion has a monotonically-increasing scale factor, a{t), the overdensitics will 
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have a scale radius, R{t), which will eventually collapse at a finite time. For instance, 
in the case where C = 0, equations (1.7.8) and (5.2.3) can be automatically integrated 
to give: 



The difference between those two solutions clearly increases, especially after the 
turnaround of the ovcrdcnsity, when that region starts to collapse. As the collapse 
proceeds, the bigger will be the difference between the badcground and the over- 
density. Variations in a between the background and the overdensities are therefore 
quite natural although they have always been ignored in studies of varying constants 
in cosmology. 

5.5 Discussion of the Results and Observational 
Constraints 

The development of matter inhomogeneities in our universe affects the cosmological 
evolution of the fine structure constant [2, 3]. Therefore, variations in a depend 
on nature of its coupling to the matter fields and the detailed large-scale structure 
formation model. Large-scale structure formation models depend in turn on the 
dark-energy equation of state. This dependence is particularly strong at low red- 
shifts, when dark energy dominates the density of the universe [161, 162]. Using the 
BSBM varying-a theory, and the simplest spherical collapse model, we have stud- 
ied the effects of the dark-energy equation of state and the coupling to the matter 
fields on the evolution of the fine structure constant. We have compared the evo- 
lution of a inside virialised overdensities, using the standard (A = 0) CDM model 
of structure formation and dark-energy modification (wCDM). It was shown that, 
independently of the model of structure formation one considers, there is always 
a spatial contrast, 5a /a, between a in an overdensity and in the background. In 
a SCDM model, 5a/ a is always a constant, independent of the virialisation red- 
shift, see Table 5.5. In the case of a wCDM model, especially at low redshifts, 
the spatial contrast depends on the time when virialisation occurs and the equation 
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of state of the dark energy. At high redshifts, when the wCDM model becomes 
asymptotically equivalent to the SCDM one, 5a/ a is a constant. At low redshifts, 
when dark energy starts to dominate, the difference between o; in a cluster and in 
the background grows. The growth rate is proportional to \Wfi\, see Tables 5.6, 5.7 
and 5.8. These differences in the behaviour of the fine structure constant, its "time 
shift density contrast" (Aa/o;) and its "spatial density contrast" {5a/a) could help 
us to distinguish among different dark-energy models of structure formation at low 
redshifts. 

Variations in a also depend on the value and sign of the coupling, C/^j of the 
scalar field responsible for variations in a, to the matter fields. A higher value of 
leads to a stronger dependence on the density contrast of the matter inhomo- 
geneities. If the value or sign of C,/uj changes in space, then spatial inhomogeneities 
in a occur. This could happen if we take into account that on small enough scales, 
baryons will dominate the dark matter density. The sign and value oi C,/uj will 
change, and variations in a will evolve differently on different scales. If there are 
no variations in the sign and value of C/u;, then the only spatial variations in a are 
the ones resulting from the dependence of the fine structure constant on the density 
contrast of the region in which one is measuring a. At first sight, one might conclude 
that the difference between etc and a;,, is only a consequence of the coupling of a to 
matter. In reality this is not so. It is always possible to obtain the same results for 
any value \C/^\ with a suitable choice of initial conditions, as can be seen from the 
results in Tables 5.1, 5.2, 5.3 and 5.4. 

The results of this chapter confirm the natural explanation, given in chapter 4, 
for why any experiment carried out on Earth [8, 165, 166], or in our local solar 
system [57], gives constraints on possible time- variation in a that are much stronger 
than the magnitude of variation that is consistent with the quasar observations [35] 
on extragalactic scales. The value of the fine structure constant on Earth, and most 
probably in our local cluster, differs from that in the background universe because 
of the different histories of these regions. It can be seen from the Tables 5.5, 5.6, 
5.7 and 5.8 that inside a virialised overdensity wc expect Aa/a ^ —10"'^ for Zy < 1, 
while in the background we have A.a/a ~ —10^^, independently of the structure 
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formation model used. The same conclusions arise independently of the absolute 
value of the coupling C,/uj, see Tables 5.1, 5.2, 5.3 and 5.4. 

The dependence of a on the matter-field perturbations is much less important 
when one is studying effects of varying a on the early universe, for example on the 
last scattering of the CMB or the course of primordial nucleosynthesis [46, 47, 48]. 
In the linear regime of the cosmological perturbations, small perturbations in a will 
decay or become constant in the radiation era [2] . At the redshift of last scattering, 
z — 1100, it was found in [3] that l^a/a < 10^^. In the background, the fine 
structure constant, will be a constant during the radiation era [107] so long as the 
kinetic term is negligible. A growth in value of a will occur only in the matter- 
dominated era [1]. Hence, the early-universe constraints, coming from the CMB 
and primordial nucleosynthesis, are comparatively weak, Aa/o; < 10~^, and are 
easily satisfied. 
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Table 5.1: ^ = — 10~^: Time and space variations in a obtained for the correspond- 
ing redshifts of virialisation, z^, for (/uj = — 10~^. Wc have assumed a ACDM 
modeh The indexes 'b' and 'c', stand for background and cluster respectively. The 
italic and bold entries correspond approximately to the level of the Oklo and (3- 
decay rate constraints, respectively. The quasar absorption spectra observations 
correspond to the values of Aa/a\i,. The initial conditions were set in order to have 
ac{z = 0) = ao- 
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Table 5.2: j = — 10~ : Time and space variations in a obtained for the correspond- 
ing redshifts of virialisation, z^^ for C^juj = — 10~^. Wc have assumed a ACDM 
modeh The indexes 'b' and 'c', stand for background and cluster respectively. The 
italic and bold entries correspond approximately to the level of the Oklo and (3- 
decay rate constraints, respectively. The quasar absorption spectra observations 
correspond to the values of Aa/a\i,. The initial conditions were set in order to have 
ac{z = 0) = ao- 
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Table 5.3: ^ = — 10~^: Time and space variations in a obtained for the correspond- 
ing redshifts of virialisation, z^, for (/uj = — 10~^. Wc have assumed a ACDM 
modeh The indexes 'b' and 'c', stand for background and cluster respectively. The 
italic and bold entries correspond approximately to the level of the Oklo and (3- 
decay rate constraints, respectively. The quasar absorption spectra observations 
correspond to the values of Aa/a\i,. The initial conditions were set in order to have 
ac{z = 0) = ao- 
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Table 5.4: j = — 10~^°: Time and space variations in a obtained for the cor- 
responding redshifts of virialisation, Zy, for (/uj = — 10~^°. Wc have assumed a 
KCDM model. The indexes 'h' and 'c', stand for background and cluster respec- 
tively. The italic and bold entries correspond approximately to the level of the Oklo 
and /3-decay rate constraints, respectively. The quasar absorption spectra observa- 
tions correspond to the values of Aa/ajfe. The initial conditions were set in order 
to have ac{z = 0) — a^. 
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Table 5.5: (A = 0) CDM: Time and space variations in a for the corresponding 
redshifts of virialisation, Zy, in the (A = 0) standard Cold Dark Matter model. The 
indexes 'b' and 'c', stand for background and cluster respectively. The italic and bold 
entries correspond approximately to the Oklo and to the /3-decay rate constraint, 
respectively. The quasar absorption spectra observations correspond to the values 
of Aa /a\b- The initial conditions were set for ^/o; = — 2 x 10~^, in order to have 
aJz = 0) = an. 
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Table 5.6: w — —0.6 CDM: Time and space variations in a for the corresponding 
redshifts of virialisation, z^, in the w — —0.6 CDM model . The indices '6' and 'c', 
stand for background and cluster respectively. The italic and bold entries correspond 
approximately to the Oklo and /9-decay rate constraints, respectively. The quasar 
absorption spectra observations correspond to the values of Aa/alft. The initial 
conditions were set for C,/uj ~ —2 x 10~^, in order to have ac{z = 0) = ao- 
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15.297 


6.362 


9.510 


4.42 



Table 5.7: w — —0.8 CDM: Time and space variations in a for the corresponding 

rcdshifts of virialisation, z^, for the w = —0.8 CDM model. The indices 'b' and 'c', 
stand for background and cluster respectively. The italic and bold entries correspond 
approximately to the Oklo and /3-decay rate constraints, respectively. The quasar 
absorption spectra observations correspond to the values of Aa/a\b. The initial 
conditions were set for (/lo — —2 x 10~^, in order to have ac{z = 0) = ckq. 



5.5 Discussion of the Results and Observational Constraints 
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14.418 


6.408 


9.469 


4.52 



Table 5.8: A CDM: Time and space variations in a for the corresponding redshifts of 
virialisation, z^^ in the KCDM model. The indices 'h' and 'd , stand for background 
and cluster respectively. The italic and bold entries correspond approximately to 
the Okie and /3-decay rate constraints, respectively. The quasar absorption spectra 
observations correspond to the values of Aa/ajb. The initial conditions were set for 
C,/uj — —2 X 10~^, in order to have ac{z = 0) = ccq. 



CHAPTER 6 



Conclusions and Future Directions 



Enquanto ndo desvendarmos este paradoxo do fim, 
por um lado o desejo de ter um telos, 
e por outro o desejo de ndo ter um fim, de ser eterno, 
continuaremos agrilhoados sobre a verdade... 

- Ernesto Mota - 

The aim of this thesis was to study the influence of matter field inhomogeneities 
on the evolution of the fine structure constant. 

In chapter 2, using a phase plane analysis, we studied the cosmological evolution 
of a time- varying fine-structure constant. We considered the case of a power-law 
a{t) — for the scale factor of the background universe. We have shown that in 
general a increases with time or asymptotes to a constant value at late times. We 
have found general asymptotic solutions for all the different and possible behaviours 
via the analysis of the critical points of the system that determines the evolution of 
a. In particular, we have also found asymptotic solutions for the dust, radiation and 
curvature-dominated FRW universe which also generalise the asymptotes found in 
[108]. These solutions correspond to late-time attractors that describes the ijj and 
a evolution in time and will be useful to study the linear regime of the cosmological 
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perturbations. 

In chapter 3, by applying the gauge-invariant formahsm of ref. [142], we have 
determined the evolution of small inhomogeneities in Sa/a in the presence of small 
adiabatic density inhomogeneities. The evolution of the perturbations of the metric 
and of the stress-energy filling the universe behave as in cosmological models with 
constant a and we can to determine the behaviour of small inhomogeneities in 5a /a 
in the gravitational fields created by the density and metric perturbations. In a 
flat radiation-dominated universe we find that inhomogeneous perturbations in a 
decay on large scales while on scales smaller than the Hubble radius they undergo 
bounded oscillations. In a flat dust-dominated universe, linear perturbation theory, 
indicate that small inhomogeneities in a will become constant on large scales at 
late times while on small scales they will increase as t^/^. In reality, the small-scale 
evolution will be made more complicated by the breakdown of the assumptions 
underlying the perturbation analysis and the development of local deviations from 
the FRW behaviour. In an accelerated phase of our universe, as is the case for 
an early inflationary epoch, or during a A- or quintessence-dominated late phase 
of evolution, we show that inhomogeneous perturbations in a will decrease on all 
scales. This result complements the earlier discovery [107], [108], that a tends 
to a constant with exponential rapidity in Friedmann universes dominated by a 
cosmological constant. Any pre-existing inhomogeneities will be frozen in but their 
scale will be exponentially increased by the de Sitter expansion. 

In chapter 4, we have shown that spatial variations in a inevitably occur because 
of the development of non-linear inhomogeneities in the universe. We used the 
spherical collapse model to study the space-time variations in a in the non-linear 
regime. Strong differences arise between the value of a inside a cluster and its 
background value and also between clusters. Variations in a depend on the matter 
density contrast of the cluster region and the redshift at which it virialised. If 
the overdense regions are still contracting and have not yet virialised, the value 
of a within them will continue to change. Variations in a will cease when the 
cluster virialises so long as it does so at moderate redshift. This leads to larger 
values of a in the overdense regions than in the cosmological background and means 
that time variations in a will turn off in virialised overdensities even though they 
continue in the background universe. The fact that local a values 'freeze in' at 
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virialisation, means we would observe no time or spatial variations in a on Earth, or 
elsewhere in our Galaxy, even though time- variations in a might still be occurring on 
extragalactic scales. For a cluster, the value of a today will be the value of a at the 
virialisation time of the cluster. We should observe significant differences in a only 
when comparing clusters which virialised at quite different redshifts. Differences 
will arise within the same bound system only if it has not reached viral equilibrium. 
Hence, variations in a using geochemical methods could easily give a value that is 
100 times smaller than is inferred from quasar spectra. 

In chapter 5, we study the dependence of the fine structure constant on the equa- 
tion of state of dark energy and on the coupling to matter, C/o; . It was shown that 
independently of the model of structure formation we consider or the coupling to 
matter, there will be always a spatial contrast, Sa/a, between a in an overdensity 
and in the background. When comparing (already) virialised overdensities, in the 
Standard Cold Dark Matter model {SCDM), spatial variations of the fine struc- 
ture constant are independent of the redshift we measure them, and are constant. 
However, that is not the case with time-variations which grow in time. A simi- 
lar behaviour is found, at high redshifts, in a ACDM or other dark-energy CDM 
models. At low redshifts, spatial variations of the fine structure constant become 
a function of the virialisation time, the density of the clusters and the equation of 
state of the dark energy component. We conclude the chapter, showing that the 
natural possibility of occurring, time or spatial inhomogeneities in C,l^i would lead 
to spatial and time variations in the fine structure constant. 

The independence of the results above, in particular the effect of the matter 
inhomogeneities on the cosmological evolution of the fine structure constant, with 
respect to the model used to drive variations in a, is an important issue. All the 
results we have obtained were only performed in the case of the BSBM theory. 
In spite of being a quite simple and elegant theory, the BSBM model but does not 
include the effects coming from grand unification theories. These may be important, 
since in many situations the field responsible for the variations of a is coupled to 
all the other matter fields. The coupling to the other matter fields might result 
on effects which are more important than the matter inhomogeneities, and could 
avoid the spatial inhomogeneity 5a/ a. Until similar studies are performed to other 
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models, which are more complex and "realistic", for instance low energy string- 
theory models, there will be always doubts about the universality of the results. 

Another important result, which requires a more detailed investigation, is our 
claim that a "freeze in" after virialisation occurs. That conclusion was taken only 
based on qualitative arguments. Hence, a quantitative analysis of the process of 
virialisation of the overdensities, and the behaviour of the fine structure constant 
inside the overdensity, is of major importance. These can be obtained performing a 
more complete numerical simulation where many more variables and features of the 
structure formation models are involved. 

If the results above are indeed confirmed in the future, and for other theoretical 
models, the consideration of the evolution of inhomogeneities, notably the one inside 
which we live, is essential if we are to make meaningful comparisons of different pieces 
of astronomical and terrestrial evidence for the constancy of a. 
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